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ABSTRACT 


Volume 1 of this report discusses the radiative and collisional 
processes which determine the recombination rate in a plasma. A gen¬ 
eral expression is derived for the recombination rate, and numerical 
values are given for a hydrogen plasma in the temperature range 250°K 
to 64,000®K and for all electron densities. These results are applied to 
an expanding hydrogen plasma and to the expanding debris from a nuclear 
explosion, and indicate that in an explosion occurring in a vacuum, sig¬ 
nificant recombination may take place at times later than a millisecond 
after the detonation. 

Volume II contains seven papers dealing with the physics of recombi¬ 
nation in plasmas. Included are the radiative transition probabilities for 
a number of atoms and ions obtained from a theory of atomic screening 
constants. 
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THE STATIONARY PLASMA* 


2.1 The Basic Equations 

In order to Investigate the non-equilibrium properties of a 
partially ionized plasma, it is necessary to consider the electronic and 
atomic processes which affect the populations of the individual atomic 
energy levels. A particular level may be populated by spontaneous and 
stimulated transitions from other discrete levels, by Impact induced 
transitions from bound states, by radiative recombination of free elec¬ 
trons and by three-body recombination processes involving free electrons. 
It may be depopulated by the inverse processes of radiative transitions 
out of the level, impact induced transitions out of the level, photo- 
ionization and impact induced ionization. 

The increase in population per second of a particular level (p) is 
given by the differences between the rates of population and depopula¬ 
tion. The equations for the populations in each of the levels will 
have the form 

^ "<P> <2-l) 

J P 

In this equation n(p) is the population of the p^^’ level. Tlie expres¬ 
sion n(J) appearing in the summation represents the rate at which 
atoms are entering the p^^ level from the level. Tlic term b^ 
represents the rate of capture of free electrons into the p^*' level 
and the expression a^^ n(p) represents the rate at which atoms are 

*Thi8 chapter is repeated in Volume II to make It available in an un¬ 
classified document. 
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leaving the level for other levels and the ionized state. The a's 
and b's include the effects of both radiative and collision processes, 
and they are functions of the density and temperature of the free elec¬ 
trons. Explicit expressions for the coefficients are developed in 
Appendix A. 

2.2 Comparison of the Radiative and Collisional Processes 

Tables giving the coefficients for the radiative and collisional 
processes are collected in Appendix A. The tables for hydrogen are 
essentially complete. For other elements the tables represent a start 
on the cookiete assembly of data needed for estimating recombination in 
complex plasmas. For the present they serve as an indication of the 
compilation and are useful in showing the pattern of variation of the 
coefficients. 

An examination of the tables shows that the spontaneous transition 
probabilities decrease rapidly with increasing excitation, in contrast 
to the collisional excitation and ionization probabilities, which in¬ 
crease rapidly. Thus, from the tables for hydrogen in Appendix A, it 
may be observed that even if the electron density is as low as 10^^ cm ^ 
and the electron temperature is as low as 1000“K, electron impact ioniza¬ 
tion of level n = 10 is more rapid than the loss by spontaneous transi- 
20 -3 

tion. For an electron density of 10 cm and an electron temperature 
of lOOC^K, electron impact ionization of the first excited level is more 
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rapid than the loss by a spontaneous transition. Electronic excita¬ 
tion and de-excitation processes are still more efficient so that a 
Boltzmann distribution of the highly excited states is to be expected 
except at low electron densities and tenq>eratures. 

2.3 Previous Calculations of Recombination Rates 
( 2 - 2 ) 

Baker and Menzel have applied equations of the form of 

Equation (2-1) to an Investigation of the relative intensities of 
spectral lines in gaseous nebulae in which the ground state of atomic 
hydrogen is ionized by stellar radiation. In these circumstances, a 
steady state prevails so that 

“ n(p) = 0 

and Equation (2-1) reduces to a set of algebraic equations, the solu¬ 
tion of which is equivalent to the inversion of a matrix. Baker and 
Menzel do not distinguish between the individual sub-levels corres¬ 
ponding to a given principal quantum number and consequently assume 
implicity that 

n(n,i) = \ n(n,i) (2-2) 

n 
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f 2-3 ^ 

Burgess^ ' has repeated the calculations taking account of the 
different sub-levels and he gives the equilibrium populations for electron 
temperatures of 10,000°K and 20,000°K. He finds that the ^-levels are 
far from being populated according to their statistical weights and that 
for an optically thin plasma, the p-levels (i-1) are much under-populated 
due to (he large probability of np - Is transitions. 

In both Investigations, all collision processes are Ignored. Al¬ 
though this Is a permissible approximation for the low-lying levels of 
hydrogen when the electron density and the electron temperature are not 
high (as Is the case In must gaseous nebulae). It will ultimately be 
Invalid as n Increases. 

Since collision processes which redistribute the sub-levels corres- 


ponding to a given n, such as 


e + H(2p) - e + H(2s) 

(2-3) 

H'*' + H(2p) - h’’ + H(2s) 

(2-4) 


are very efficient, the effect of Including collision processes will 
tend to ensure Equation (2-2). The more refined treatment of Burgess Is 
therefore actually less accurate than that of Baker and Menzel when the 
electron temperature Is high and when the principal quantum number n 
Is large. 

Glovanelll^ ^ has Included electron collision processes In his 
treatment, which again assumes statistical equilibrium. In his first 
paper In which he assumes Equation (2-2,))ic shows that Saha's ionization 
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formula applies to the populations of the high-lying levels, provided the 
electron teiiq>erature is enq>loyed. This occurs because the most efficient 
depopulating mechanism is electron impact ionization and the most effi¬ 
cient populating mechanism is the inverse process of three-body recom¬ 
bination and these processes are very rapid. The same situation must 

apply even in non-equilibrium cases. 

(2-5) 

In a second paper Giovanelli discusses the populations of the 

low-lying levels. The study contains an inconsistency, which may be 
serious, in that while he distinguishes between the individual sub- 
levels, he ignores redistribution processes such as Equation (2-3). 

His results are therefore invalid at high electron densities and tempera¬ 
tures. We note in passing that his results are also invalid at low elec¬ 
tron densities and temperatures in that he ignores the decay of the 
metastable 2s state by the two-quantum process 

H(2s) - H(ls) + hvj + hV2 (2-5) 

Discussions of the hydrogen plasma in non-equilibrium situations 

(2-6) (2-7) 

have been given by D'Angelo' ' and by Bates and Kingston' . In 

both investigations, heavy particle collisions are ignored and it is 

assumed that the plasma is optically thin. Both investigations assume 

Equation (2-2) and ignore processes such as Equation (2-3). This is 

equivalent to supposing the Equation (2-2) is very fast. 

D'Angelo obtains an approximate solution of Equation (2-1) by attemp¬ 
ting to follow through in detail the reaction paths of any individual 
electron, making the assuiqption that recombination is effective only for 
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those paths which terminate in the ground state. He Ignores the effect 
of electron impact excitation and de-excitation. He gives results for 

three electron temperatures, lOOO^K, 3000®K and 10,000°K and for electron 
12 13 -3 

densities in the range from 10 to 10 cm only. His Investigation 
reveals the importance of three-body recombination processes (as did the 
earlier work of Giovanelli) but his method of solution is laborious to 
apply and the accuracy of the numerical values of the effective recombina¬ 
tion coefficients is uncertain. 

The method employed by Bates and Kingston is simple and more accurate. 
They argue that the rates at which atoms in a particular excited state 
are produced or destroyed are very much greater than the rate at which 
the very small number density of atoms in that state changes. It follows 
that for an optically thin plasma except for the ground state 

^ n(p) - 0 (p . 2, 3, ...) (2-6) 

so that recombination is effective only into the ground state. Then 
writing 

S"(‘) 

(where n(l) signifies the ground state) 

defines the effective recombination coefficient which will 

depend upon the electron density as well as the electron temperature. 

Their siethod reduces the set of coupled differential equations, Equation 
(2-1), to a set of coupled algebraic equations. The number of levels re¬ 
quiring Individual examination is restricted since the populations 
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of the highly excited levels are given by the Saha formula, the level 
at which this becomes true being provided by the calculations. 

2.4 The Number of Discrete Levels 

In none of the investigations we have described has the question of 
the number of levels that should be Included in the equations been exam¬ 
ined. It is clearly not infinite because ultimately the discrete levels 
will overlap the continuum. 

/ 2-1 \ 

Ivanov-Kholodnyi, Nikol'skii and Gulyaev'' ' have discussed the 
lowering of the continuum due to pre-ionization. The treatment is not 
entirely satisfactory but should suffice for our purposes. They find 
that the number n^ of levels of different principal quantum number which 
remain discrete is given by 

21.65 - £n(n ) 

£n(n^ + 1) -g-^ (2-8) 

Values of n^ have been conq>uted from Equation (2-8) for a wide range 
of electron densities and they are presented in Figure 2-1. The values 
are quite large and it seems probable that Stark broadening of the lines 
will cause a greater effective depression of the levels. The Stark 
broadening has been discussed by Griem,Kolb and Shen^ and it may 

be necessary to carry out detailed computations based on their theory. 
However, the fact that the highly-excited levels are in essentially 
thermal equilibrium and play only a small role in determining the effect¬ 
ive recombination suggests chat the solution of the equations is insensitive 
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("u) 31VJ.S QNnoa 1S3H9IH 

Figure 2-1. Highect Bound State of an Electron i.n Hydrogen 
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ELECTRON DENSITY (cm'^) 



to the nuniber of levels Included. The Stark broadening will be Inqpor- 
tant, however. In determining the opacity, which In turn effects the 
recombination. 

2.5 Computation of Recoiriblnatlon Rate In Optically Thin Hydrogen 

(2-7) 

Adopting the method of Bates and Kingston , we first change 
the variables In Equation (2-1) from n(p) to d(p) by Introducing the 
relationship 

n(p) - n^ n(X‘^) 0 )^ exp(l^/kT^) e(p) (2-9) 

(where m is bmss of electron, k Is Boltzmann's constant, h Planck's 
constant and 1^ the Ionization potential of this p^^ level.) 

The new variables d(p) are the ratios of populations In the p^^ levels 
to the populations derived from the Saha relationship. At equilibrium 
d(p) Is unity. 


The Introduction of (2-9) into Equation (2-1) gives an equation of 
the form* 


^e(p) .jp «(1) 

jVp 


b 

P 


a 

PP 


e(p) 


P = 1, 2, 3, ... 


(2-10) 


Assuming with Bates and Kingston that 

MM ,,0 p = 2, 3, 4,... (2-11) 


*the a's do not have same meanings In Equation (2-10) as In Equation 
( 2 - 1 ) 
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we obtain 


° “X ®jp ‘’p “ ®pp 

jTp 

p = 2, 3, 4,... (2-12) 

with the equation for the ground state where p ^ 1 treated as a special 
case: 


^ e(l) = y Sjj e(J) + bj - a^j 0(1) (2-13) 

jVp 

If the series of equations in (2-12) are truncated at some level k, the 
Equations (2-12) are in the form of k-1 linear equations with k unknowns, 
the unknowns being taken as the 0's. By inversion of the matrix, a 
solution may be obtained for all 0's in terms of one unknown, 0(1). Be¬ 
cause of the linearity of the equations, the solution will have the form 

0(p) - C^ + Cj 0(1) 

p = 2, 3... k (2-14) 

The right side of Equation (2-13) can then be evaluated also in terms of 
the single unknown 0(1) and taking account of Equations (2-9) and (2-7), 
one suy obtain the effective recombination coefficient in the form: 

■ *. - *1 

Values of a^ and a^ in Equation (2-15) have been obtained for an 
optically thin hydrogen plaatu by oiachine methods for electron tempera¬ 
tures ranging from 2S0'’K to 64,000'*K and for all electron densities. 
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11 


indices give the power of 10 by which the entries must be multiplied. 
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Indices give the power of 10 by which the entries must be multiplied. 



levels up to n = 10 being Included in the confutation. Preliminary 

values of a and a, ha« been provided by Bates and are given In Tables 
o X 

2-la and 2-lb, respectively. 

There are certain Irregularities In the values as given In the 
tables which suggest there might be an error In parts of the tabulation. 
Until a further check of the computations has been coofleted the values 
In the table should be regarded as uncertain by a factor of 2. 

2.6 Dy[Ci»^gj,op.,o| Cojapjitatlpfl flRd SSSHUg 

In order to Illustrate the significance of the differential 
Equations (2-10), the algebraic Equations (2-12), and the solution as 
given In Tables 2-la and 2-lb, It will be worthwhile to consider the 
changes which take place In the populations of the various excited levels 
of a hydrogen plasma. Suppose that for an optically thin hydrogen plasma, 
the temperature and number density of the free electrons are maintained 
constant by some arbitrary means. Under these conditions, the rate of 
change of the populations of all the bound states are given by the 
differential Equations (2-1), and for given Initial conditions which 
assign a value to the population of each of the bound states at the out¬ 
set, a solution of the Equations (2-1) will give the populations of each 
of the states for all subsequent time. The general character of this 
solution Is Illustrated schematically in Figure 2-2. 

In the figure, the populations of the several states are described by 
the variable 0 , which Is the ratio of the number In a state to the 
number In that state under conditions of thermodynamic equilibrium. 

Since we are supposing that the number density and temperature 
of the free electrons reiiulns constant. 
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RATIO OF POPULATION TO EQUILIBRIUM POPULATION, 


0IDA40- 1310 



TIME -- 


Figure 2-2. Scheaatic Dlagraii of Population Changes During 
RecoMblnaclon 
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6(p) will be proportional to n(p). In Figure 2-2, it is supposed that 
the ground state is unpopulated at the outset, but the initial popula¬ 
tions of the various excited states are taken at values which have no 
particular significance. As indicated in the figure, at the outset there 
will be a very rapid change in the populations of the upper states. In 
particular, the value of 6 for the upper states will change nuch store 
rapidly than that of the ground state. After a short time, the rapid 
changes in 6 of the upper states will decline, and a pseudo-steady state 
will occur such that, if viewed on a short time scale, all the populations 
remain constant. There will, however, be a non-zero rate of increase in 
the ground state, which will siake itself felt over a longer period of 
tisie. This pseudo-steady state condition is representated in Figure 2-2 
as occurring Just before the first tisie break, and it is this condition 
which is described by the algebraic Equations, 2-12 and by their solu¬ 
tion obtained taking the recombination coefficient or - a^. At an appre¬ 
ciably later period of tisie, as represented after the first time break 
in Figure 2-2, the Increase in 6 for the ground state will be appreciable, 
and there will be a corresponding change in the pseudo-steady state values 
of the populations of the various higher states. This situation is re¬ 
presented in the middle time-period of Figure 2-2, and is also described 
by the solution to the algebraic Equations, 2-12, but now since 6(1) is 
not equal to zero, it is necessary to take the effective recombination 
rate, as o - a^ - aj0(l). 

The populations of the various states continue to change as recomb¬ 
ination proceeds, and the right hand section of Figure 2-2 indicates 
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the general nature of these changes on a much different time scale than 
the portions of the figures at the left and the center. The populations 
will eventually approach asymptotic values at which they are constant, 
and at which time the net recombination ceases. This condition, which 
Is indicated at the extreme right of the figure represents the steady 
state under conditions of radiation loss. It can be seen that for this 
condition, the population of the ground state exceeds that which It would 
have for thermpdynamlc equilibrium. The population of the ground state 
can be obtained readily from the Tables 2-la and 2-lb by setting a ■ 

a^ - aj 0(1) - 0. 

In simplifying the differential Equations 2-10, to the corresponding 
algebraic Equations,. 2-12, one gives up any possibility of obtaining 
Information about the Initial variation shown at the extreme left of 
Figure 2-2. Except for this Initial response, which Is unimportant for 
the applications considered here, the effective recombination rate 
and the populations of all the states can be obtained throughout the 
period Illustrated In the second and third time blocks In the Figure 
from the algebraic equations. 

It should be observed that in the proceedure followed here, one 
supposes that the temperature and number density of the free electrons 
and the conditions on light absorption as described by the X'P^i'enBters 
are given conditions of the problem. To do otherwise would make the 
equations (either the differential Equations 2-10 or the algebraic 
Equations 2-12) non-llncar and would vastly Increase the difficulties of 
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obtaining a solution. The solution obtained here can nonetheless be 
applied to important problems in which the temperature, density and 
light absorption in the plasma change with time because the pseudo¬ 
steady state to which Equation (2-12) applies is achieved in times 
short compared to those for significant variation in the plasma properties. 
This matter will be discussed further in Section 3 of Volume I. 

2.7 PUfma O^h^r. Than Hydyogen 

The general plasma has several ionized components each of which 
can exist in several states of ionization. Nevertheless, the general 
plasma may be treated in a manner similar to that employed for the hy¬ 
drogen plassia, but with certain modifications. 

Just as for the hydrogen plasma, it is necessary to collect to¬ 
gether all the radiative and collislonal transition probabilities. Un¬ 
fortunately, no systematic tabulation of these probabilities has been 
given even for ground states of systems and the collection of the data 
is a major task. We have made a beginning in assembling the data and 
the tabulated value are included in Appendix A. 

The basic equations are similar to (2-1) and their solution pre¬ 
sents no problem provided heavy particle collisions can be ignored. 

With this proviso, there is a set of equations of the form (2-1) for 
each component of the plasma and these sets of equations are uncoupled. 

The quasi-steady state method may then be applied to each component 
separately. 
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It does, in fact, appear reasonable to neglect heavy particle 
collisions except for charge transfer processes involving the ground 
states. This exception arises because, although charge transfer into 
the ground state of a system may be inefficient, it does couple compo¬ 
nents together and it has a permanent recombining effect. It does not 
introduce any difficulties into solving the equations since the ground 
state equation receives special consideration in any event. 
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APPENDIX A 

A. Dalgarno, A. Naqvi, E. Ashley and H. E. Stubbs 
AT(»11C PROCESSES IN A RECOMBINING PLASMA 

The time dependence of the population of each state of a recombining 
plasma can be obtained from a consideration of the several processes by 
which the atoms enter and leave a given state. Both radiative and 
colllsional processes must be considered. 

A.1 Radiative Transitions Between Bound States 

The rate of depopulation of the state p by spontaneous transitions 
to lower states q, may be written 

n(p)^ A(p,q) (A-1) 

q<p 

where A(p,q) is the probability of a spontaneous transition from the 
state p to the state q measured in sec 

If the plasma is not optically thin, stimulated emissions and absorp¬ 
tions occur. For the rate of upward transitions which depopulate the 
state p, we write 

^ ^ X(p,r) A(r,p) (A-2) 

’’ r>P 

where is the statistical weight of the state p and X(p,r) is a factor 
which takes account of Che opacity. In particular for a plasma in ther¬ 
modynamic equilibrium with the radiation field 



X(r,p) 


<exp (hv 



(A-3) 


where energy of the r - p transition. For an optically thin 

plasma,Xvanishes. The rate of stimulated downward transitions can be 
written 


n(p) 


X(p,q) A(p,q) 


(A-4) 


The p state can also be populated by bound-bound radiative transi¬ 
tions and we have for the rates of population by the three processes 


^ n(r) A(r,p) 
r>p 


t>p 


X(r,p) A(r,p) 


(A-5) 


(A-6) 


“o Z ' (q.P) A(p,q). (A-7) 

<i<p q 


The probability of a spontaneous transition from a level p to a 
lower level q may be expressed in the form 


-.43 

A(p,q) - — y ' " S(p,q) 
3hc u 
P 


(A-8) 


where v is the frequency of the emitted radiation, u is the statistical 

P 

weight of the p^^ level and S(p,q) is the transition strength, defined 
according to 

S(p.l) - .^ |< pjT, > I (A-9) 
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where'P and are the configurations of the initial and final states 
^ P q 

and r is the position vector of the j " electron. In the case of Russell- 
Saunders coupling, S(p,q) reduces to the form 

S(p.q) -/(L) /(M) (A-10) 

where /(L) is the relative strength of a multiplet within a transition 
array, ^(M) is the relative strength of a line within a multiplet and 
a is a transition Integral involving radial wave functions only. Tables 
of C^lL) and (M) are available for most multiplets (Goldberg^^ 

White and Eliason^^ Russell^^ , Allen^^ and if necessary, they 
can easily be extended following Rohrlich^^ The transition integral 

0 presents greater difficulties. It R(ni|r) is the initial radial wave 
function of the active electron and R(n'i'|r) is the final wave function, 

2 

(A-11) 

where is the lower of the azimuthal quantum numbers and Z* which 
characterize the active electrons. 

For hydrogen, values of a^(nX,n'i!^ have been computed by Green, 

Rush and Chandler^^ and the corresponding values of A(n£, n'jJ*) for 
n and n* less than or equal to 8 are presented in Table A-1. Of interest 
also are the total transition probabilities A(n,n') obtained by summing 
k{n£, n*Z') over all possible values of i and i' consistent with n and 
n', assuming that the sub-levels are distributed according to their 
statistical weights. The total transition probabilities so obtained are 
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presented in Table A-2. There are ocasional differences from previous 
tabulations of Menzel and Pekeris^^ 

For more complex systems, approximations must be made. For many 
cases and especially for the excited states of light elements, the 
Coulond) approximation of Bates and Damgaard^^ should be adequate. In 
this approximation the active electron is regarded as moving in the 
field of an effective nuclear charge being determined from the observed 
ionization potentials for the states. Bates and Damgaard have presented 
2 

tables from which a may be determined for transitions involving small 
values of i! and / and we have used them to calculate the transition 
probabilities of the more Important lines of Li I, Li II, C I, C II, 

C III, N I, N II, N III, 0 I, 0 II and 0 III. The results, illustrative 
of the information which may be obtained by the Bates and Damgaard 
methods, are presented in Table A*3. Also, Table A-4 lists the results 
of calculations by this method of the transition probabilities of 
resonance lines of a number of atoms of interest. 

For heavier elements it may not be possible to apply the Coulomb 
approximation since the ionization potentials are not always known and 
a less empirical approximation scheme is required A theory using 
screening constants, which is described in the following appendices, 
is capable of dealing with these cases. 
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DROGEN IN 


8 9 10 ' 11 12 13 


» 

i) 

0 

0 

») 


> 

) 

) 

;) 

) 


1.232 (5) 
4.674 (4) 
2.299 (4) 
1.286 (4) 
7.799 (3) 
3.007 (3) 
3.357 (3) 
2.329 (3) 
1.663 (3) 
1.215 (3) 
9.065 (2) 
6.882 (2) 


7.147 (4) 
2.811 (4) 
1.426 (4) 
8.187 (3) 
3.007 (3) 
3.323 (3) 
2.267 (3) 
1.597 (3) 
1.156 (3) 
8.530 (2) 
6.449 (2) 


4.374 (4) 
1.773 (4) 
9.226 (3) 
5.414 (3) 
3.422 (3) 
2.278 (3) 
1.577 (3) 
1.126 (3) 
8.248 (2) 
6.170 (2) 


2.797 (4) 
1.163 (4) 
6.182 <3) 
3.876 (3) 
2.376 (3) 
1.605 (3) 
1.126 (3) 
8.136 (2) 
6.023 (2) 


1.856 (4) 
7.880 (3) 
4.268 (3) 
2.595 (3) 
1.692 (3) 
1.158 (3) 
8.219 (2) 
6.002 ( 2 ) 


1.270 (4) 
5.493 (3) 
3.024 (3) 
1.865 (3) 
1.232 (3) 
8.527 (2) 
6.116 (2) 


14 15 16 17 18 19 


8.928 (3) - 

3.923 (3) 6.425 (3) 

2.191 (3) 2.863 (3) 4.717 (3) - 

1.368 (3) 1.619 (3) 2.128 (3) 3.528 (3) 

9.139 (2) 1.022 (3) 1.217 (3) 1.609 (3) 2.679 (3) 

6.394 (2) 6.899 (2) 7.763 (2) 9.291 (2) 1.234 (3) 2.066 (3) 
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TABLE A-3 


BOUND-BOUND TRANSITION PROBABILITIES A(p,q) FOR THE IMPORTANT 
LINES OF Lil, LlII, Cl, ClI, CIII, NI, Nil, NIII, 01, 011 AND OIII 

___ __ Lll Lines _ 


Transition 

^2 

•^1 

X(X) 

2 

CT 

A(10^sec 

2p^P°-28^S 

1/2 

1/2 

6709.8 

5.424 

3.63 


3/2 

1/2 

6709.6 

5.424 

3.63 

3p^P°-28^S 

1/2 

1/2 

3233.6 

0.01641 

0.098* 


3/2 

1/2 

3233.6 

0.01641 

0.098* 

4p^r°-28^S 

1/2 

1/2 

2742.0 

0.01392 

0.136* 


3/2 

1/2 

2742.0 

0.01392 

0.136* 

5p^P°-28^S 

1/2 

1/2 

2563.0 

3.90x10’^ 

0.0468* 


3/2 

3/2 

2563.0 

3.90x10"^ 

0.0468* 


Lill Lines 


Transition 

x(X) 


2 

a 

AdO^sec’S 

l82p^P°-l8^ ^S 

'.v9,3 

6 

3.72.:10‘^ 

... ‘1 

1 .90* 

l83p^P°-U^ ^S 

178.0 

6 

9.647x10'^ 

691 

l84p^P°-l8^ ^S 

!_ __ 

171.6 

6 

3.219x10'^ 

257 


Cl Lines 


2p(^P°)3s^p"-2p^ ^P 

1657 

18 

- 

0.01397 

1.24* 

2s2p^ ^D°-28^2p^ ^P 

1561 

30 

00.9563 

101 

2s2p^ ^P°-28^2p^ ^P 

1329 

18 

0.01284 

2.21* 

2p(^P°)48^P°2p^ ^P ■ 

1280 

18 

0.001488 

0.287* 


^below 907. maximum value 


2R 







cm Lines 


28(^S)2p^P°-28^ 

977 

6 

0.5723 

248 

28(^S)3p^P°-28^ 

386 

6 

0.05579 

392 

28(^S)4p^P°-28^ 

310 

6 

0.006529 

88.5 


NI Lines 


2s^2p^(^P)38S- 





2,V *S° 

1200 

12 

0.007382 

0.860* 

2p^(^P)48S-2p^ ^S° 

964 

12 

0.0007751 

0.175* 

2p^(^P)3dS-2p^4S° 

953 

120 

0.04008 

93.7 

2s2p^ S-2s^2p^ ^S® 

1135 

1135 

0.5631 


2p^(^P)58S-2p^ V 

911 

12 

0.0003190 

0.0853* 


*Below 9071 maximum value 
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TABLE A-3 (continued) 
Nil Lines 
















TABLE A-3 (continued) 


on Lines (C.A.) 


Transition 

uX) 

afw 

2 

ff 

AdO^sec S 

2s2p^ S-2s^2p^ V 

834 


0.2670 


2s^2p^(^P)3sS- 





2s22p^ V 

539 

12 

0.006351 

8.2* 

2p^(^P)3dS-2p^ ^S° 

430 

120 

0.005050 

128 j 

2p^(^P)3dS-2p^ ^S° 

430 


0,005356 



01II Line s 


2s2p^ ^D°-2s^2p^ 

3p 

834 

30 

0.1417 

"" "1 

98,4* 

2s2p^ ^P°-2s^2p^ 

3p 

703 

18 

0 1538 

179 

2s2p^ ^S°.2s^2p^ 

'P 

508 

12 

0,1643 

1014 

2s^2p(^P°)3s^P®- 
2s V 'P 


374 

18 

0 006293 

48.6* 

2s^2p(^P°)3dV- 
2s22p2 3p 


306 

135 

0,005453 J 

346 


'‘'^Belov 90% maximum value 
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TABLE A-4 

TRANSITION PROBABILITIES FOR RESONANCE LINES 


Element 

Transition 
(Upper - Lower) 

''z -'^1 

X(X) 

cr^ 

AdO^sec" 

Rb 

I 

4p^(^S)5p 

4p^(^S)58 ^S 

1/2 - 1/2 

7947.60 

8.82 

35.5 

Sr 

I 

58(^S)5p ^P°- 
58^ Is 

1 - 0 

4607.331 

4.72 

195 

Sr 

II 

4p^(^S)5p V- 1 
4p^(^S)58 ^S 

1/2 - 1/2 

4077.714 

4.92 

147 

Zr 

I 

4d^58(aS)5p ^F°- 

4 - 4 

6121.49 

3.19 

52.5 



4d^58^ ^F. 

2 - 2 

6134.58 


49.6 

Zr 

II 

4d^(a^F)5p V- 5 
4d^(a^F)58 S 

1/2 - 4 1/2 

3391.96 

3.59 

186 

Nb 

I 

4d^58(a^F)5p ^F°- 
4d^(a^D)58 S 

5 1/2-4 1/2 

4058.933 

3.39 

410 

Mo 

I 

4d^(a^S)5p ^P°- 

4 - 3 

3798 359 

3.35 

123 



4d^(a^S)58 ^S 

3-3 

3864.115 


117 




2 - 3 

3902.968 


114 

Tc 

I 

5d^58(a^S)5p S°- 

3 1/2-2 1/2 

4298.2 

3.15 

160 



Sd^S.^ 

2 1/2-2 1/2 

4263.4 


164 




1 1/2-2 1/2 

4239.4 


167 

Ru 

I 

4d^(aS)5p V- 
4d^(aS)58 ^P 

Entire 

Multiplet 

3799.124 

3.27 

120 
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TABLE A-4 (continued) 


Transition ' *^1 « ? fi -l 

Element (Upper - Lower) X(a) a AHO^sec ) 


Rh 1 

4d®(^F)5p V- 

3 1/2-4 1/2 

3692,357 

3,21 

115 


4d®(^F)58 

2 1/2-3 1/2 

3657 987 


1082 

Ag I 

o 

CO 

1 1/2- 1/2 

3280 682 

2.95 

169 


4d^°(^S)58 

1/2 - 1/2 

3382.890 


154 

Cd I 

4<i'°5«O^S)5|> *P“ 
4d“>5.2 ‘S 

- 1 - 0 

2288 7 

1 99 

671 

Tn I 

58^(^S)68 ^S- 

1/2- 1 1/2 

4511 310 

0 890 

39 2 


58^(^S)5p V 

1/2 - 1/2 

4101 764 


26,1 

Sn I 

58^5p(^Pjl)68 

’-1-2 

3175 046 

0 202 

10 6 


2 2 3^ 

58^5p^ ^P 

0 - 1 

3034 120 


29.2 

Cs I 

5p^(^S)6p ^P®- 

1 1/2- 1/2 

8521 10 

9 94 

32.5 


5p^(.^S)68 

1/2 - 1/2 

8943 50 


28 1 

Ba 1 

68(^S)6p ^P®- 
68^ 

1 - 0 

5535 484 

5 78 

138 

Ba 11 

5p^(^S)6p ^P®- 

1 1/2- 1/2 

4554 033 

5 96 

127 


5p^(4)68 

1/2 - 1/2 

4934,086 


100 

U I 

5d68(a^D)6p V- 

2 1/2-2 1/2 

7270.11 

5,10 

50.0 


5d68^ 
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TABIf A-4 (continutd) 


Element 

Tranaltlon 
(Upper - Lower) 

J - J 

2 ‘'l 

Hh 

0^ 

AdO^aac’^ 

Ih 1 

4d®(^P)5p V- 

3 1/2-4 1/2 

3692.3S7 

3.21 

115 


4d®(^P)5i S 

2 1/2-3 1/2 

3657.987 


1082 

Ag I 

4d^°(^S)5p ^P®- 

1 1/2- 1/2 

3280 682 

2.95 

169 


4d^°(^S)5a 

1/2 - 1/2 

3382.890 


154 

Cd I 

4d^°5aO^S)5p ^P® 

- 1-0 

2288.7 

1.99 

671 


4d'°5.2 ^ 





In I 

5a^(4)6e V 

1/2- 1 1/2 

4511.310 

0.890 

39 2 


5a^(*S)5p V 

1/2 - 1/2 

4101.764 


26.1 

Sn I 

5a^5p(^pji)6a ^P®- 1 - 2 

3175.046 

0 202 

10.6 


5a^5p^ ^P 

0 - 1 

3034.120 


29.2 

Ca I 

5p®(4)6p V. 

1 1/2- 1/2 

8521.10 

9.94 

32.5 


5p^(^S)6e 

1/2 - 1/2 

8943.50 


28.1 

Ba I 

6a(^S)6p ^P®- 

1 - 0 

5535.484 

5.78 

138 


6a^ 





Ba 11 

5p^(^S)6p ^P^- 

1 1/2- 1/2 

4554.033 

5.96 

127 


5p®(^S)6a 

1/2 - 1/2 

4934.086 


100 

U 1 

5d68(a^D)6p V- 

2 1/2-2 1/2 

7270.11 

5.10 

50.0 






A. 2 lUdlatlw Bound-Fr«< Tramittona 


In order to calculate the coefficient a(p) for radiative capture of 
electrons into state p 

+ e ^ X(p) + hv 

we must make some assumption about the electron velocity distribution. 
Since the electrons exchange energy readily in elastic collisions, they 
may be described by a Maximilian distribution characterised by an elec¬ 
tron temperature The recombination coefficient is accordingly 

“<p|V*7(7y^<"*V 

X J %i?) (hv)^ exp(-hv/kT^) [l +X(v)] d(hv) 

** (A-12) 

vhere 1^ is the ionisation potential of the p^^ level, hv is energy of 
the radiation photon and Q^(p) is the photoionisation cross section 
corresponding to the process 

X(p) + hv - + e (A-13) 

The rate of population by radiative capture is 

a(p|T,) n(X*) (A-14) 

vhere n^ and n(X^) are respectively the electron and X^ ion number den¬ 
sities. 


35 


The rate of depopulation by photoionization may be written 


n(p) P(p) 

where ^ 

p(p) - f 

The quantal description of bound-free transitions Is entirely similar 

to that of bound-bound transitions except that the radial wave function 

R(n'f'|r) in (A-11) becomes R(€'C'|r), describing an electron of positive 

energy e* moving In the field of the Ion. Thus, the cross section for 

the transition of an electron from the bound state (ni!) to the continuum 

state of energy e' Is 2 

Ancea 

Q^(ni!;€’) - 

£'•&-! 

where is the Ionization potential, C^, are simple algebraic quantities 
which depend upon the Initial and final atomic configurations and 
g(nX;e'/') Is a radial transition Integral 

g(ni;€’i') ■ r^R(n.2jr) R(e'i|r)dr (A-17) 

o 

Although exact analytical expressions are available for the matrix 

elements describing the photoionization of the excited states of hydrogen 
(A-9) 

(Gordon they are complicated and have not yet been given a numerical 

form except for the lowest lying levels. Giovanelli^^”^®^ employs an 
approxlisate formula due to Gauntfor the photoionization cross 


tjjl 



I 
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which in numerical form is 




2 . 82 x 10 ^^ 
5 3 


(cm‘) 


(A-19) 


P V 

This choice has also been made by Elwert^^~^^^ and Ivanov-Kholodnyi 
et but it is of questionable accuracy. Recently, Burgess 

has presented tables useful for low velocities of the ejected electron, 
from which more accurate photoionizetion cross sections may be obtained 
for a hydrogen plasma. His calculations distinguish between the various 
sub-levels and apply to all levels with p £ 12. 


Using his calculations, the radiative recombination coefficients may 
be calculated according to (A-12) and Burgess has given values of a(nA|T^) 
for electron temperatures of 10,000°K and 20,000^K. He do not wish to 
distinguish between the individual sub-levels (except possibly in the 
latest stages of the plassw decay). For a hydrogenlc plasma of nuclear 
charge Z, after suitably averaging over the sub-levels, we obtain the 
formula(Seaton 

a(p,Z|T,) = DZ Xp^^2 Sp(X) (A-20) 

where 

D . 5 ^ ( 5 )' “f ' *0 ' 5 W-21) 

(A.22) 

and 

Sp(X) * (Xp) + X’^^^ (Xp) + x‘^^^ (Xp) (A-23) 
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with 

(x) exp (x) Ei(x) 


(A-2A ) 


(x) 0.1728 x^^^ (u V l)‘^^^(u - (A-25) 


(X) - 0.0496 x'^^ r (u I t I u + Ue’^^^clu (A-26) 

o 


nM 

p" 


r 


■ dv (A-27) 


Seatoi/^'^^^ has published lablus ol xS^‘^^(x), xS^^^.^x) and xS^“\x), iroro 
which the radiative recombination coeiticienls a(p,2lT^) may be calculated 
readily. Values obtained lur hydrogen (Z 1) in this way have been provided 
by Bates and ate shown in Table A-^. Seaton also presents tables similar 
to those for the xS^'^^x) Irom wliicli the total recombination coefficient 
to all levels 


‘-I I %■ 

p'sl 

may be calculated directly. 


(A-28) 


For more complex systems the only method generally applicable is an 
extension of the Coulomb approxinution. Tills yields formulas identical to 
(A-16) and (A-17) except that n is replaced by an effective (,uanlum number 
n* - z/V 1^^> X being the residual charge on the atom when the electron 
has been ejected. 
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The indicles give the powers of 10 by which the entries oust be niltiplied. 














According to Burgess and Seaton 

gCn’^r ;e V ) . ^ ) cos n [n* + (i(e') + X(n*i;e'X')] (A-29) 

C(n*J) 


G(n*i;€'?') 




Itt 'n* ^jL£' 


(A-31) 


Iffc'n* 

n(e') is the continuum phase shift obtained by extrapolation of the quantum 
defects (i(I .) of the energy levels of the bound states. 


?(n*, i) 


2 -N.ii) 

be 


(A-32) 




• (n* > f + 2) 


tabulated by Burgess and Seaton for 1 S n* £ 12 
and /, /' <3. For values of n* < i ♦ 2, additional tables are presented 
which permit the calculation of g(n*i;€'i') for transitions from specific 
bound states to the continuum. 


As an example of the application of Burgess and Seaton's 
formulas and tables to the treatment of complex atomic systems, they have 
been used to calculate the photoionization cross sections and radiative 
recombination coefficients for some of the observed states of aluminum. 
The ground state configuration of Al I is 3s^3p and the observed excited 
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2 

states, with the exception of one 3s3p configuration, arc configurations 
of the form 3s^nf where 3 < n < 11 and 0 £ f < 3 (cf. C. Moore^. The 
latter configurations consist of a single (nf) electron outside a closed 
shell for which the C^, assume the particularly simple form (selection 
rules on the azimuthal quantum number allow only transitions where 
f -- f ± l)i 


jg-t-l ^ 

‘^i+l " 2i+l ’ ^jC-I " 2f+l 


(A-34) 


Table A-6 gives the Burgess and Seaton coefficients for the allowed 
transitions In A11. These coefficients were used to calculate the func¬ 
tional parameters individual transitions 

The results, together with the Ionization potentials and quantum defects 
of the bound states listed in Table A-7. The ^(n‘' ,£), also 

given in Table A-7, were determined from the curves of Figure A-1 which 
shows plots of the continuous phase shifts for each(nX) series. 


For small positive electron energies e', the phase shift may be taken 
as a linear function of e'. The extrapolated quantum defects determined 
from the curves of Figure A-1, lead to the following approximate expres¬ 
sions for u(e'): 


H(€') = 1.755 


(ns series) 


H(£') = 1.264 

- 0.065e' 

(np series) 

(A-35) 

H(£') = 0.985 

+ 3.3c’ 

(nd series) 


H(€') - 0.044 

+ 0.26c' 

(nf series) 
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fhotoionizatlon eroit ••ctions ware calculated for all but the 3a3p^ 
and n£ levels, for ejected electron energies c* lying in the interval 
(Kydbergs) 0 ^ e'— 0.5 for the nd levels, 0 5. c'l. 0.3 for the ns levels 
and 0 1 e'— 0.2 for the np eeries. For electron energies in excess of 
the maxiDum values indicated, it was found that the aj^roxisMte foraulaa 
led to incorrect results. This siay in part be due to the aseusiption of 
a linear dependence for the however, it should be noted that the 

approximations Inherent in the derivation of (A>29) limit the validity 
of the calculations to small ejedted electron energies. 

The photoionization cross sections obtained are listed in Table A-8. 

2 2 

The 3s 3d and 3s 4d cross sections were calculated under the assumption 
that n* S.*2 since explicit expressions are not tabulated by Burgess and 
Seatonfor nd “* e'p snd nd “* e'f transitions when n* < >4+2. The 
resulting cross sections are therefore of lesser accuracy than the remain¬ 
ing cross sections shown. 

The radiative recoad>ination coefficient amy be calculated from Bq. (A-12) 
once the photoionization croas sections arc known as functions of the ejected 
electron energy. Thus, in terms of electron energy for a thin plasma 

• 7 ( 

Instead of evaluating the integral of (A-36), the recombination coefficient 
can be determined approximately by use of the ei^ression 

a(n/|T^) » V (lU) (A-37) 


4f 




b* ■ultiplisd. e* is Im Ifydbsrgs, 



















































where (ni) the crose section for radiative recombination is given by 
detailed balancing 

(I i+€)^ 

Qg (nf) - (2i+l) —“- (ni) (A-38) 

me € 

and V Is the average electron speed 

/.« 

V «■”> 

Substituting (A-39) and (AOS) in (A-37), we find 

- ^ (VO^Q„(n.) 

. , Q («/) 

- 2.046x10® (2i+l) (lp+€) ~i/2 


The radiative recoarivination coefficients calculated from (A-40) are 
given in Table A-9 for electron teoperatures of 6190®K and 24,760*K. 



RADUTZVE REGOMBZHATION COEFFICIENTS 
FOR A1 I a(ni|T Xn cm^/gec) 


- 6190®K Tg - 24,760®K 


s 

s 

P ' 

d 1 

;_^_i 

8 

P 

d 

3 


2.11(-13) 

6.48(-14) 



8.04(-15) 

1.94(-14) 

4 

9.80(-15) 

3.90(-14) 

|2.36(-14) 


1.99(-16) 

1.73(-15) 

1.53(-14) 

5 

1.86(-16) 

1.32(-14) 

1.04(-16) 


5.34(-17) 

7.68(-16)i 

2.56(-15) 

6 

7.28(-16) 

5.62(-15) 

9.39(-16) 


9.24(-17) 

2.68(-16) 

6.80(-16) 

7 

4.94(*16) 

5.06(-15) 

1.33(-15) 


1.01(-16) 

2.96(-16) 

2.21(-16) 

8 

1.74(-17) 


1.28(-14) 


B.35(-17) 


1.31(-16) 

9 

8.35(-18) 


1.10(-16) 


7.63(-17) 


1.03(-16) 

10 

2.22(-18) 


8.66(-16) 


B.04(-17) 


8.46(-16) 

11 



7.33(-16) j 




7.52(-16) 


The brackets give the power of 10 by which the 
entries must '/O multiplied. 






A.3 Electron Collision Induced Tr«n»ltlon» B«tween Bound 


Let (Ptq) be Che cross section for the excltstlon to the q stete 
by Che impect of sn electron of energy e with en atom in the p state. 

The rate of depopulation of tht p^*' state by electron impact is 


r > p 


(A-41) 


where the collision probability C(p,r) depends upon Che electron temp¬ 
erature and is obtained by integrating Che product of the electron 
velocity and the cross section over the Maxwellian distribution accord¬ 
ing to 

C(p,r) « expC-c/kT^) dc (A-42) 


The probability of de-excitation to a lower level is similarly 


, n(p)^C(p.q) 

q 


(A-43) 


3 




The direct and Inverse collision probebllities are related through the 
principle of detailed balancing: 

U) f 1 - 

C(p,q) - 'if' ‘xp -^JC(q,p) (A-44) 


The p state may also be populated by electron intact induced 
excitation and de-excitation and the appropriate rates are 


.. I 

q < p 

.. I 


n(q) c(q,p) 


n(r) C(r,p) 


(A-45) 


(A-46) 


r > p 

There are various simple ways of computing the probid>ilities of 
excitation and de-excitation transitions induced by electron-atom 
collisions (c.f. Allen^A-"^ none of which is entirely satsifactory. 
Those quantum mechanical treatamnts in which some reliability may be 
placed (e.g. Mott and MasseyBates et al^^”^^\ Milford and co- 
workers^^'McCrea and Me Klrgan^^’^^\ and Me Carrol^^*^^^ have 
been based on the Born approximation, and therefore, are valid only 
for high electron impact energies. 


SI 



A Blnpllfled formila daacrlblng tha crosa-aectlon for •xeitatlon froH a 
q level to a p level Induced by the form (c.f. Allen^^*^^) 

QJq.p) - 1.44 X lo“xa^ (ch)^f(q.p) jC 7 ^ - 7 ) (c«^ (A-47) 

<!*? 

where e denotes the lispactlng electron energy, and ^ is the ex¬ 
citation energy, both expressed in Rydbergs (c^ ^ ■ c^- . f<q,p) represent^ 

the oscillator strength for absorption, and is related to the spontaneous 
transition probability A(P>q) by 

3 u 

f(q.p) - —fVj (4-48) 

8 x C q 

Equations (A-47) and (A-48) auty be used to calculate collisional transition 
probabilities for the atomic systems for which the ipontaneous transition 
probabilities are given in Tables A-1, A-2, and A-3. If Q^(p,q) denotes 
the corresponding collisional de-excitation cross section, then (A-47) ' 
and (A-48), with the aid of the following expression (c.f. Allen^^*^^): 

Vp * Vq ^€^‘**P^ 

«p*(p.q) - «^f(q,p) (A-50)' 


yields 


Q^(P.q) - 5.78 X 10“ xe 


f(q 




(A-51) 


Ixpoeasions sisdlar to (A-47) and (A-Sl) with valuee of A(p,q ) obtained 
by the formulas of Msnsel and Pekeris^^*^^ for hydrogen were employed 
by Ciovanelli^^"^®^ and Ivanov-Kholodnyi et 


S2 


Perhaps the most accurate method of computing colllslonal transition 
probabilities is the classical method of Gryzlnski^^’^^^, which represents 
an improvement over the original formulation of Thomson^^'^^\ in that 
it takes account of the fact that the atomic electrons are moving. The 
cross-sections so desired fall off too rapidly at high electron energies, 
as i ^ instead of the correct (c ^ In e), but their accuracy j^obebly 
suffices, the effect of excitation and de-excitation not being a pri¬ 
mary one. 


Considering in sosm detail the classical transfer of momentum during 
a coulomb collision between an incident electron of energy e and an atomic elec¬ 
tron having binding energy Grysinski derived the following approximate ex- 
prcti.sior for the differential cross-section aU)t) corresponding to a 
loss in energy of the incident electron of magnitude Ac; 


o(<ife) « 




3 ’" 


.fSft < € - 6 


1 - 




(A-52) 

The total cross section for a collision involving an energy loss ^ subject 

to the restriction € < c < e is then given by 
p — ~ q 

P 

Q-<€ ,€ ) • fa(ae)d(aie) • Q(€ ) - Q(c ) (A-53) 

€ p q J p q 

where Q(e^) the cross section for a collision with electron loss greater 

than Is defined by Equation (A-^5). Thus idien ^ - e , as in the case 
P P 

for Erogenic orbits» 


1/2 
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2 « < « 
P 


^1 “ (cm^) 12 C 2 e 


(A-54) 


Thtt collision trsnsition priAsbility is given by 

\U2 

M S’"*- 1 O <, , , 

“ q’ p 




(A-55) 


Inserting *p) leads to the following collision prob¬ 

abilities for hydrogen, where ■ l/p^ end - 1/q^ 

C(P.,) - ^ ^ A •■' y| 

e q i qp 


(A-56) 


(A-57) 






3 [2 -(g» y(l » 4g)] (y - l)^^^(g -f 1) 


6^'* (g^y)^'^ 


^C-3g ♦ y(3 ♦ 4g)] y'^^ 


7g + y)*'* 


y < 1 ♦ g 
(A-5« 


y > 1 + g 
(A-59) 


C(p,q) has been calculated for hydrogen by Bates using equation 
(A-56). The results are given in Table (A-10) fcnr p & q less than 
or aqual to 10. 


54 



1 

1 

I 

I 

! 


i 

"« 

S 
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tout A-10 (coatlntwd) 



1 h« teMkats give the powers of 10 by iriiich the entries mst be sniltiplled. 

























A-10 (continusd) 



Hm tar«ck«ts give th* powers of 10 by which the entries wust be siultiplied. 
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For non-hydrogenie syetems having known energy levels the more 
general expressions of (A-52) and (A-53) must be used for the calcula¬ 
tion of collision probabilities. 
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.4 Electron ColllBlon Induced Bound-Free Trantltloni 


Let q(p,e) d« be the cross section for the formation of electrons 
with energy In the range c, c + de by collision of an electron having 
energy E with an atom In the p state. The total Ionisation cross section 

Is 

E-Ip 

Qe(p) - i J qgCP.O (A-60) 

and the probability of Ionization by collision Is 


"(to) J exp(-E/kTg) (A-61) 

The rate of depopulation of the (p) state due to electron Impact lonlza> 
tlon Is 

n(p) % D(P) (A-62) 


The Inverse process Is three-body recombination: 


+ e + e ^X(p) + e. (A-63) 

The efficiency of the three-body process may be written In the form 
m • 

FCp) J J «2 exp(-«j/kT^) expC-fj/kT^) 




(A-64) 


ifhere y(p|<^C 2 ) probability that an atom Is found In the (p) state 

as a result of the collision with a positive Ion and two electrons with 
energies and C 2 respectively. The probability is related to the in¬ 
verse process of electron li^iact ionisation by detailed balancing. Thus 
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16 Tma). 

<Ip + «1 + *2> " “5- * 1*2 Y(pI*i. «2^ 

P 

Where ou^ Is the statistical weight of the continuum electron and 

B - Ip + Cj + Cj . (A-66) 

The rate of population of level (p) by the three-body recombination 
process (Equation A-63) Is then 

n^ n<X+) P(p) (A-67) 

Alternatively, the probability of three-body recombination F(p) may 
be obtained by a slispler method noted by Ivanov-Kholodnyl, et al ^^ 

They remark that In thermodynamic equilibrium, the number of particles 
leaving a given quantum level via any specific process Is equal to the 
number of particles entering this state via the Inverse process. It 
follows that 

P(p) - cp(p) D(p) (A-68) 
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vrtiere 


, Anvk. T \-3/2 

cp(p) . expdp/kT^) (A-69) 


As in the case of excitation, the cross sections decrease too rapidly 
as the energy of Impact Increases. Since Ionization and Its Inverse 
are the controlling processes, the error may be serious and we have 
therefore given some attention to the problem of how we may obtain Im¬ 
proved accuracy, whilst still retaining the universal applicability and 
the simplicity of the Gryzlnskl method. A new universal method that 
has resulted Is described below. 


Again for reasons similar to those mentioned above for bound- 
bound collision probabilities, the most accurate method for computing 
electron Impact Ionization probabilities, which Is of universal ap¬ 
plicability Is that of Gryzlnskl^^"^^\ Giovanelli^^”^®^ and 
Elwert^^"^^^ used the classical Thomson fonoula which can be expected 
to give somewhat less accurate results than that of Gryslnskl. 
Ivanov-Kholodnyi et al^^~^^^ employed approximate expressions simi¬ 
lar to those for the collisional excitation cross-sections. 


»:3 


The classical Thomson cross-section (cf. Thomsonor 
Fowlerfor the ionisation of an atom having n^ valence elec¬ 
trons by an electron with energy €, with the corresponding ejection 
of an electron from a p-state into an energy interval C + d(t is 


Q, (P.O di = - di 


(A-70) 


The Thomson total ionisation cross-section is then 


4 ns a. At 

% <P) = —V-®- / ; 

‘ ‘ ''o «+ ey 




(A-71) 


a^ is the radius of the first Bohn orbit in hydrogen. 
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Ai pointed out in the previous section, the Grysinskl formulss 

represent sn laprovestent over the ThoiMon derivation in that accouht 

is taken of the motion of the bound atomic electrons. According to 

Grysinskl the total ionisation cross-section for an electron-electron 

collision in vhich the impacting electron suffers an energy loss AE 

greater than the threshold ionisation energy s is given by 

’’ 

Q(. ) . I <?(4D d (IS) (A-73 

S 

where is the maximum permissible transfer of energy consistent 

with conservation of momentum and energy, and a (AE) is the differ¬ 
ential cross-section described by Equation (A-52). (The expression 
for Q(c) given in (A-54) %fas obtained from the integration indicated 
by (A-72)). Thus the cross-section for ionisation from the p level 
is Just (A-S4). The ionisation transition pr(d>ability for a Haxwel- 
lian electron energy distribution is given by 

d. U-73 


Making the substitutions x - e/s^ we find from (A-S4) and (A-73) that 



(cm^/sec) X > 2 




X < 2 


In the ceee of hydrogen, ■ 1/p^ fo thut 



Values of D(p| T^) in hydrogen (provided by Bates) are given in Table 
A-11 for p less chan or equal to 10. For the more general ease, Equa¬ 
tion (A-74) must be used to compute D(nl | T^). 

As in the case of excitation, Che cross-sections decrease too 
rapidly as the energy of impact increases. Since ionisation and its 
inverse are the controlling processes at the higher energy levels, 
the error may be serious and we have therefore given some attention 
to Che prcAlem of how to obtain improved accuracy, %rhilsc still re¬ 
taining the universal applicability and the simplicity of the Grysin- 
ski method. 


Bewricing Equation (A-Jl) for the niomson ionisation cross-section 
in the form 


<«) 


,2 1 


K*-?) 


«id introducing the diswnsionless parameter \ • €/l 


(A-?6) 


66 

















(A-47) becomes 


I Oi'^n . 21A 1 > 

—T*-- T' It ) 


(A-77) 


^e rlght-hend side of (A-77) is e universel function of independ¬ 
ent of the perticuler ionising process. 


It is known from the asymptotic form of the Born epproximetion 
(cf. Mott end Masseythat for large values of X, behaves 
asymptotically according to 

Qj - X*^ In X (A-78) 

so that (A-77) is incorrect in the limit of high impact energies. 

The modification of (A-71) by Grysinskialso has this defect. 


Ue shall attempt to improve (A-77) by introducing a correcting 
factor f (X) such that 




f(X) 


(A-79) 


4* a ‘ g(X) 


and determine f (X) by eppeel to experimental data. 

In figure (A-2), all the available experimental data^^”^^ 
on the electron impact single ionisation of atoms has been 
used to determine l\\ •• a function of X, adopting for n^ the 
mniber of electrons in the valence shell. If f (X) exists such that 
(A-79)is true, this procedure would produce a single curve. Instead, 
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there results distinct curves which are similar in shape but have an 
appreciable spread in magnitude of the croas-scctions. This spread 
representa the posaible error that would be a consequence of the use 
of the mean of the curves of Figure (A-2) as a universal function. 

lie may reduce this factor by taking partial account of the elec¬ 
tron structure of the atomic systems, it being an oversimplification 
to adopt for n almply the nusd>er of valence electrons n^. The signifi¬ 
cant quantity is actually the equivalent number of dispersion elec¬ 
trons for transitions into the contlnum. 

two alternative ways. We may 
the same magnitude at large 
that the maximum value of 

(A-80) 

(A-81) 

The latter procedure is the simpler and we adopt it. It leads to the 
curves shown in Figure (A-3), all of the normalised experimental data 
for i<mitation of neutral atoms lies within the shaded area. No 
curve now departs from the mean of all the curves by a factor of more 
than about lOX (cf ahaded area of Figure A-3) and we shall select the 
mean as a universal curve. The normalised curve for ionisation of 


To determine n, we may proceed in 
choose n to ensure all the curves have 
energies or we may observe from (A-71) 
is given by 


and determine n from 
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Figure A-3. Uaivereal Curve for Electron Impact Ionization 
Cross Sections 



of He^ Ion falls Just outside of this area and Is also shovm In Figure 
A-3. 

The corresponding values of n are given In Table (A-12a). We see 
that there Is a pattern to n which allows us to predict It with some 
reliability at least for the lighter eloaents. Thus, since n appears to 
be an Increasing function of atomic number, we can assign Integer values 
of n to appropriate ranges of the atomic number as shown In Table (A>12b) 

From the universal curve of Figure (A-3), we can proceed to the de¬ 
termination of a universal Ionization probability curve as a function of 
electron temperature T^. Thus from Equation (A-S8) we find 

D(Te) = (A- 82 ) 

Using Equation (A-77) 

D(Te) = (kT^)“^^" 4« na^^ jT* Xg(X)e‘^^*'^e jX (A-83) 

where T^ = T^/1 Is the electron temperature measured In units of the 
Ionization potential. Finally 

0(j^)T^3/2 ^ 4ji na^^y'*Xg(X)e"^^’‘^e dX (A-84) 

The right hand side of which Is a function of the electron temperature 
T^ (Kily and may be computed once for all atomic systems for any given 
value of T^, using the universal curve g(X) given In Pisure (A-3). 
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lASLI A-12a 

EQUIVALBMT MUHBBR OF DISIltSXON KUCTBOMS, 
n KTBRNiMn) Wi BQOAXlOli (A-ll) 








DT 

The function . Is shown as a function of T' In Figure (A-4). 

n e 

The Ionisation probability for any element for any temperature may be 
read off the graph, the appropriate abscissa being obtained by dividing 
T^ by the Ionisation potential I. 


1 

I 

! 
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Figure A>4. Ifciiveraal XonlMtloii ProlMiblllty Curve 


I 
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A.5 Heavy Particle Colltitong 


Consider now Cransltlona produced by collisions with heavy particles 
which we label Y^. Let (p.r) be the cross section for the excitation 
to the r^^ state by the iaipact of a heavy particle with an atom 
X in the p^^ state, the energy of relative sntion being e. Hie rate of 
depopulation of the p^^ state by heavy particle, impact is 


"(P) Z Z 

I r>p 


where the collision probability Pj.(p,r) is obtained by multiplying the 
cross section by the relative velocity of the collision and integrating 
over ion velocity distribution f(v) to give 

Pj(P.r) ■J V (p,r) f(v) dv (A-86) 

If f(v) is the Maxwellian velocity distribution corresponding to temp- 


( (p,r) € exp(-€/kT) de 


where p is the reduced mass. 


The probability of de-excitation to a lower level is similarly 


n(p) P^(p.q) U-Sa) 

i q<p 
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The rates of population are 


^ n(Y^) ^ n(q) (q.p) (A-89) 

1 q<p 


^ n(Y^) ^ n(r) (r.p) (A-90) 

I i>p 


The ionisation and recombination involving collision with ions may 
be formulated in a manner similar to that for electron collisions. The 
rate of depopulation by ion impact induced ionisation may be written 


n(p) { mx *) Q(p) + Z "(Yj) Qj(P)) (A-91) 

CD 

Q(P) - —(kT)*^^^ J Qgip) exp (-B/2kT) dE 


(XM,)’ 


(A-92) 


Q^(p) is Che total ionisation cross section of the state (p) by ions of 
energy B. 

The rate of population by three>body recon^ination 

X"*” + Yj + e -► X(p) Yj (A-93) 

may be written 

nCx”*”) n(Yj) t(p) (A-94) 

where R(p) is an expression exactly similar to Equation (A-64) 


77 


since for allowed trensltione the croes sect lone for the Impact of 
an electron and a singly charged positive Ion of the s«ne velocity are 
of much the same size, excitation and Ionization by Ion liqpact Is usually 
very Inefficient compared to electron Impact. For exaa^le, three*body 
recombination 

X'^ + X"*" + e -► + X (A-95) 


will be negligible compared to 

x‘‘‘ + e + e-*X-fe 


(A>96) 


E;:ceptions may occur when the energy defect of the process Is 
small as In 


h"*" H(n,i) h'*’ + H(n,/') 


(A-97) 


and for transitions between pairs of highly excited states. The effects 
of Ion Impacts Is to strengthen the argument that a BoltssMnn equili¬ 
brium prevails for the excited states. The position regarding charge 
transfer reactions Is analogous, charge transfer being rapid at low 
velocities only when the energy defect Is low. 


We merely wish to provde a few nuatbers In order to show quantitatively 
that, apart from re-distributlng processaa, (and transitions between 
pairs of highly excited states), electron collisions arc swre efficient 
than heavy particle collisions. 
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(A-34) 

Seaton' has calculated the rate coefficient for the re-distribution 

process 

+ H{2s) h'*' + H(2p) (A-98) 

o -4 3 

at 10,000 K, and he obtains the very large value of about 5 x 10 cm 

sec Re-distribution may also occur by charge transfer 

h"*" + H(2s)-m H(2p) + h'*' (A-99) 

but Boyd and Oalgarno^^ have shown that (A-99) is much less rapid 
than (A-98). 

Excitation and ionisation proceed much more slowly and may be 
neglected. Thus, the maximum values of the cross sections for the 
direct processes 

h'*’ + H(ls) -► H* + H(2s) (A-100) 

h'*’ + H(ls) - + H(2p) (A-101) 

are respectively about 5 x 10 cm^ at about 20 kev impact energy, 
(Bates^^ and 1 x 10 cm^ at about 30 kev (Bates^^ and 
the cross sections are decreasing rapidly with decreasing energy. 
Similarly, it may be seen from the calculations of Bates and 
(A'SS A-39) 

Dalgarno' ' that excitation accosipanying charge transfer such 

as 

h'*’ + H(l8) H(2p) (A-102) 
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Is extrsmely inefficient for temperatures below several kev. 

Thus, for a hydrogen plasma or other singly Ionized plasma of a 
single component, processes due to Ion Impact may be Ignored. When more 
than one component exists In a plasma, however, it may be necessary to 
consider heavy particle processes, for charge transfer processes are 
unique In coupling together the two components. This fact may be seen 
from Equation (A-106), by noting that If the heavy body coefficients 
?!, Q^> and are set equal to zero, all the terms containing n(Y^), 
the 1 .umber density of other Ions vanish. Account must be taken not 
only of ordinary charge transfer 

- X"*" + (A-103) 

but also of radiative charge transfer 

X**^ + y'*' •* x"*" + y"*^ + hv (A-104) 

which may be more repld than the previous reaction. 

An Investigation of the possible rates of these cherge transfer 
reactions, one of the participants of which Is an Ionic system in Its 
ground state. Is necessary, and as an example of the latter process 

He”*^ + H He'*' -► He'*' + hv (A-105) 

may be treated with high accuracy. 

Special features also arise In the case of a single cosqwnent plasma 
In which the cosqtonent may exist In more than one state of Ionization. 


A.6 Th« Basic Eauattoni 


The increase in population per second of a particular level (p) Is 
given by the difference between the rates of population and depopulation 
that we have written in the earlier sections. Thus 

[ Z ^^**’‘*^ X(p.q)} + ;7 5 X(p,r)A(r,p)] 

^ q<p ^ 

"e C(p,r) + Z C(p,q) + D(p)J 

+ Zn(Y^) ^lp^(p.r) + EPjCp.q) +Qi(p)j 

+ 3(p)| n(p) 

•^|A(r,p) [l + X (r,p)] + n^ C(r,p) + ^ n(Y^)P^(r,p^ n(r) 
", C(qjp) n(Y^)P^(q,p)^ n(q) 

cKp) XCn”**) ♦ n^ x"*" y^Cp) 

+ He ^ "i<^t> *i^P> (A-106) 

We arrive at a set of coupled, differential equations of the form 

+ « „(p) . y , n ♦ b (A-107) 
dt pp Z- pj J p ' 

J ^ ? 

where the a's and b's can be identified by coa^arison with the 
preceding equation. 
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APPENDIX I 
A. Naqvl 


ATOiXC tlUNfSniON PROBABILITIES 
B.l INTRODUCTION 

Consider two snorgy levels, 1 end 2, of sny ston, end let their 
nsss^ies, E^ end B 2 , be sueh that E 2 > E^. Traasitiens frea 2 to 1 result 
in eaission, snd the reverse transitions, froe 1 to 2, are, aaong other 
processes, caused by the absorption of radiation, of frequency V 2 j^ given 
by 

^ (B-1) 

The asMunt of energy (in erg ce*^ sec'^) esdLtted in a spectral line 
is given by 


(B.2) 

where 02 is the number of etoms per cm^ in level 2. A(2,l) end B(2,l) 
are, respectively, the probabilities (in sec'^) of spontaneous transitions 
and the transitions stUsilatad by the presence of an external radiation 
field. I 2 j^ is the intensity of this external radiation field at the 
frequency V 2 j^. 

Similarly the amouat of energy (in ergs cm’^ sec*^) absorbed in a 
spectral line is given by 
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fii ~ BO>Ol2\ ^ 34 \ 


(B-3) 


«h«r« B(l,2) is the probability of tho upward transition, from laval 1 to 
Isval 2, with tha absorption of radiation of fraquancy V 2 |^* 

Tha thrao transition probabilitias ara ralatad to aaeh othar by 
tha folUnring ralatioas dua to Binstain. 

(B-4) 




(B-5) 


whara and §2 ara tha statistical waights of tha two lavals. It is 
tharafora sufficiant to datarmina oaa of tha thraa transition probabilitias, 
usually A(2,l). 

Tha thaory for calculating atoadc transition probiAilltias f<nr 

alactric dlpola radiation is wall known (Condon and niortloy^'^\ 

(B>2) 

SlatMr' Savaral raviawa havo boon written in racont yaws also 

(6arstang^^^^ Maqri^*^^). Tha probability A<2,1) of a apontanoous 
radiati^ tranaition is giwwi by 


M 



(B-6) 


g 2 It the ttatlttlcAl wtlght of lovtl 2, X it the wavelength of the 
spectral llna enltted, S(2,l) is known as the line strength and for LS 
eouplii^ Mn bo wrlCtaa aa 

SU,').SC<,i).TCv-)r(«)(r’- 

%ihere J*<« and f<M) are factors which arise as a result of integratim 
over the angular part of the wave function. ^(L) la known aa the rela¬ 
tive atrength of a line within tha aultiplet to which it belonga and 

ia known aa the relative aultiplet strength. Both of these quanti¬ 
ties can be calculated in a clooed fona for any given line and any given 
aultiplet. f(L) has been tebuUted by Mbite and BlUaon^^*^^ and by 
8usaell^®‘‘\ f(M) has been tebuUted by OeUberg^"^ These 

tables have also been reproduced by Allen^*®\ Becently Rehrlich^®”^^^ 
has derived general eiq^ressUns, based on Racah'a aethod, froa which 
aultiplet strengths for traneitione not covered by Goldberg can be obtained. 

The electric-dipole truMitioa integral a® is given by 

(r\n,t,, 0-8) 

idMre the P'a are r tiaee the radial wave foMtioaa, and is the greater 
of the two asiauthal quaatua nuabsrs and The caleuUtion of the 
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•btolut* transition probabilitiss thoraforo d^anda upon tha ealculatim 
2 

of tha transition intatral <x » ohidi in turn dopnda an appropriata 
ehoiea of tha radial nova functions. 

B.2 THE SADUL HAVE FUMCT1(»IS 

Tha radial nova functions, P(ni;r), ara tha solutions of tha radial 
part of tha non«ralativistic Schroadingar aquation, 

whara E is tha aigan valua of tha Haailltonian for a given stationary stats 
and V(r) is tha central field in uhich tha electron novas. Hota that 
atonic units ara used throughout this report. 

Tha iaportonca of • knowledge of tha radial wave functions cannot 
be over anphasisad. Tha radial wave function, in a nutiAell, gives us 
an inaight into all other physical properties of tha aten. A knowledge of 
tha radial wave functions is requirod not only in the calculation of tha 
transition probabilities but also photoisoiaation cross sections, elastic 
scattering cross-sections, collision cross-sections for excitation of 
atsns, x-ray scatterii^ and nunorous nolecular and solid state properties, 
Moag others. 

Tha nost csnnsnly used nathod is tha self-consistent field nathod, 
developed by Hartree and extended by Fock to include exchange effects (for 
a recent review, see HartreeThese are obtained by nunwical 
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solution of s systsm of simultaneous, inhomogeneous, non>linear, integro 
differential equations, the so-called Fock equations. Efforts have also 
been issde to include the effects due to configuration mixing and correla¬ 
tion of electrons. (Sec Hartree^^^^^), The procedure is very laborious, 
and has to be performed separately for each configuration of each atom. 
Since numerical methods are used the results are given in a tabular form, 
and cannot be used for further analytical calculations. It is, therefore, 
not possible to develop analytical methods of great generality, but every 
problem must be solved separately. Thus although the theory was developed 
in all essential details over 30 years ago, a list of atosis and their con¬ 
figurations for which the self-cmsistent field radial wave functions arc 
available is pitifully small. In particular it is clear that if one wants 
to calculate atomic parameters, such as transition probabilities, for a 
large number of atosis one must develop methods of considerable generality 
which can, at one time, be applied to a large mmbmt of eases. 

On account of their obvious usefulness, many atteapts have been 
made to obtain analytic wave functions (for exas^le Zener^^^^^i 
EckertSlater^®“^^^ Morse, Young and Haurwits^®*^*^ Tubis^®'^®^ 
Norse and Yilmas^*"^^^; Breene^*"*^* Sechs and 

Weiss^^^^). Sosm of these were merely attempts to obtain analytical 
expressions of which the parameters were determined by cosparison with 
the tabular wave functions of the self-consistent field method. 

Since the advent of quantios theory, it has been very tempting to 
use the exact results obtained for hydrogen, for the more complex atoms 
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also, by th« •ubatltutlon of althar an "affactiva" principal quantum 

lUMibar n* for n, or a "aeraonad" nuelaar eharga (■ Z • s^) for tha 

Of 01 

actual nuclMT cbarga Z. Hora a^ danotaa tha acraaning of M alactron a 
by all tha aloottoma to tho atom. Tha oerooatog paramotmra hovo 

boon datarminad usually by a aami'anvirical nathod (Pauling^**^^^; 

Pauling and Shaman^^^^h, although a varUtional traatnant haa boon 
applied to aoma ainpla caaas. Tha intaraat in tha usa of hydrogonic imva 
functions with a acraanad nuelaar charge haa baan ravivad through a nathod 
for calculating tha acraaning paramatara davalopad by Laytor 

1.3 THE HYDIOGElllC HAVE FUNCTIOMS 


Tha acraanad hydrogonic radial wava functions, briefly rafarrad to 
at tha liydrogonic wavo fimetiont**, ara tho solutions of Equation 1-9 with 
V(r) raplaeod tha Geulomh potontial duo to a scraonod nuelaar charge 

E^ (■ Z • o^g), nhora s^j is a scraoning paraaotar and Z is tho atonic 
nunbor. Thoy are given by 


PCntZ^I ;r)= 



.fUf 


wlMro n* donotoa (n-i-1), tho nunhor of nodoa of tho fimetton omi |P^ la 
a eanfluont hyporgoomatrle tonotian daftood balaw. l a e a u aa of tho w^ 
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occurs In Equation B-10, wo shall find it convmiant to substituta in 
tha abova aquation. 




(B-ll) 


idiarc If' aquals 2Z/n and (> Zs^/n) is a raduead seramii^ poranatar. 


Wa now hava 




(B-12) 


1/2 

fl. r_M2L_l 


. 

iF, ^ , 2 ('♦•a 4 hi r 


(B-13) 


(B-14) 


idiora 





(B-15) 
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It li to b« noted thet different for different ni end that 

the redial wave functions (Equatione B*10 or B-12) for the sane n but dlf~ 
ferent i are not strictly orthogonal. It Is possible to nake corrections 
due to the lack of orthogonality of the wave functions, but they are not 
considered in this report. 

For subsequent use In calculating the Slater F integrals we write 
below an expression for the square of the radial wave function (Equations 
B-12 and B-U). 

[pWy^ ■,r)] = r 

^ ^ b 

(B-U) 

where 

^ 

B.4 RESUME OF LAYZER'S SCREENING THEORY 

Layser^^''^^^ has shown that the eigen values 
the correct quadratic and linear Z*dependence If the 
wave fimctions satisfy the condition that, 

P (nlyjr)= p” {nly.O+OMyiO 


of energy will possess 
one-electron radial 


(B-18) 
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where the euperecript H denotes e hydrogenic wove function end 0(n/ ^:r) 
denotes e term of order one In end 

■,r) 

so thet, et leeet, for lerge “y, 0(ni y it) term cen be neglected. Hots 
thet our notetion differs slightly from thet used by Leyser. 



Neglecting 0(ni y it) terms, end epplylng the verletlonel principle 
Leyser hes shown thet the screening peremeter n^, which gives the screen^ 
ing of the electron a by ell the other electrons In the etom, Is given by 


4 JV, 

- 


(B-20) 


where denotes the metrix elements of the Couloeib Interectlon energy 
end denotes the number of electron in the eubehell a. He cen use 
severel different epproKlmetlons for the Coulomb Interectlon, for exesple 
we mey or mey not Include cemf Iguretlon sdLxlng end we my either consider 
the everege energy of e given conflguretlon or we my edd to it the con> 
trlbutlons which come from different ^ectroscoplc terms (cherecterised 
by quentum nuiSbers S end L), belonging to the seme conflguretlon. 


In the work reported here, we heve not considered conflguretlon 
mixing, but we consider the energy oi Individual tMms es well es the 
events energy of the conflguretlon. When required, the screening con- 
stents which ere dependent on SL cen be lebelled ej^llcltly es 
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whereas those based mi the average energy will be labelled as it 


ni* 

The term energies as functions of the screening paresieters ere 

given by 


WC-n) <»-«> 

where the susnetion ie over ell the subshells in the atom. This is e very 
simple fonaule for calculating the term energies from the screening perem- 
eters. 


B.3 THE COULCMB IMTERACTKHtS OF ELECTRONS 

The most conseonly uaed epproKimotion for complex atosis (see, for 
example, Condon end Shortley^^^^) is the sopereble wove function approx- 
imetion in which each elactron is asauawd to be in a stationary state in 
the field of the nucleus and all the other electrons. The wave function 
of an N>electron atoaM is exprassed in the form of a determinant, amde up 
of the wave functions of individual electrons, aa follows, 



VI, (o') u^C«') . 

vi/a') 

\ifih •••• 


V4,Cof') Mxifli') - 



(B-22) 






I 

I 


where u is the weve function of the eingle electron problem, e^ denotes 
the eet of quentum numbers (n^ ^ m^ m^) end subscripts of u, from 1 to 
N, denote the electron to which this weve function belongs. The deter¬ 
minants! form Insures the antisymmetry property of the weve functions. 

For e meny-electron weve function of the above form, Slater 
has shown that the matrix elements of the mutual Coulomb Interactions of 
the electrons, matrix elements of the corresponding 

two electron problem. 


Ignoring configuration mixing, the average electrostatic energy 
of a configuration nay be written as 


V,. 




(B-23) 


where 




and the two electrM Interaction (a,^} Is given by 


(B-25) 


where the prlsM over the sumsHttion sign Indicates that k does not assume 
all integral values, but for all even integral values and for either 
all svm or all odd integral values, denotes the lesser of and i^. 
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Th« CMffleltnti And 8^(i^ i^) aria* fron tha angular 

and tha spin part of tha wava functions. Thaaa can ba caleulatad by the 
vail known laatrlx nathods. For details of these calculatlona tha reader 
la referred to Condon and Shortlay^^'^^ and to the more powerful nathods, 
due to Racah^^'^^\ based on tensor operators. These coefficients are 
listed In Tables B«1 end B-2. 

k k 

F end G ere Integrals over the radial part of the wave functions, 
also known as Slater Integrals. They are defined by 


. Pf^r’.r,) P(Hh Jr, <*•«> 


with 




CB-27) 



TABLE B-l 


COEFFICIENTS OF SUTER INTEGRALS FOR EQUIVALENT ELECTRONS 


_flL_^ 

I 

1 

I 

I 


_flL-flL 

0 

-2/25 

-2/63 

-4/195 


£J 

_SL-flL 

0 

0 

-2/63 

-2/143 


_SLJL 

0 

0 

0 

-100/5577 


COEFFICIENTS OF SLATER INTEGRALS FOR NON-EQUIVALENT ELECTRONS 
(Fbr « tabulated coa^inatlon (i , i.) all non-vaniahlng coaf- 
flelanta art llatad) ^ 



0 0 I -1/2 

0 110 
0 2 10 

0 3 10 

1 1 1 - 1/6 

12 10 
13 10 

2 2 1 - 1/10 

2 3 10 

3 3 1 


1 

2 

3 

4 

5 

6 

S 

I 

S 

8 

8 

8 

0 

0 

0 

0 

0 

0 

-1/6 

0 

0 

0 

0 

0 

0 

-1/10 

0 

0 

0 

0 

0 

0 

-1/14 

0 

0 

0 

0 

-1/15 

0 

0 

0 

0 

-1/15 

0 

-3/70 

0 

0 

0 

0 

-3/70 

0 

-2/63 

0 

0 

0 

-1/35 

0 

-1/35 

0 

0 

-3/70 

0 

-2/105 

0 

-5/231 

0 

0 

-2/105 

0 

-1/77 

0 

-50/3003 
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and dtnoca respectively the smeller end the gfester of the 
two radius vectors, r^^ and r 2 > of the two electrons. 

The electrostatic energy for the different terms belonging any 
configuration can now be written as 



where for some of the sln^ler cases Is as follows: 

(1) All shells completely filled, or there Is a single incompletely 
filled shell with either one electron or one hole. In these cases the 
configuration consists of a single term and, 

\4 

(ii) There la a single incompletely filled shell (designated by 
n^ /j^) with m electrons 

Note that the screening parameters are functions of the partlals 
of the Slater Integrals and with respect to the appropriate screened 
nuclear change and and G^, which we evaluate using the hydrogenlc 
wave function (Equation :'-12), arc thamaelves functions of the screened 
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nuclear charges, and therefore of the screening paraaeter it^'s. It Is 
tharefore readily seen that the rlglit«hand side of Equation 2-20 Is a non¬ 
linear function of all the screening parameters, it^, .«dien 

N Is now the nuinber of distinct subshells of electrons characterised by 
n and i quantum nuBd>ar8. 


I 

I 
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B 6 THE SIATER INTEGRALS 


From Equations (B>26). (B*27) and fB>28) it Is easy to see that 

Ic k 

F and G are synsietrlc in the two electrons. However, on Mcount of 
the occurrence of t^ and In this equation, It Is necessary to break 
up the Integral Into two parts, which for the F Integral we write as 
follows: 


4 [wi ^ I > Y I y i.) 



Fj^ = j f dr,. 

= r*'[ P(ney;r)]^ 


(3-32) 


(B-33) 


Note that F^*^ and 


are not symnstrlcal In the two electrons. 


It Is easy to prove that 
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(B-34) 


For the G Integral, theae two parta turn out to be equal and we have 

Je> Jo ^ (B-35) 


The calculatlona of the Slater F and G Integrala, uaing hydro- 
genic wave functions, is extrenely laborious, for electron shells whose 
wave functions possess more than one or two nodes. The author is not 
aware of any published work where these have been calculated as functions 
of the screened nuclear charges of the two electrons involved, nor are 
any general formulae for calculating them available in literature. We 
have, therefore, developed general formulae which are suitable for hand 
computations as well as computations using digital computers. 




<8-36) 


-H k+t K -K-l 


f, % A. Aa 


(b-37) 
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y. in,42Wi4K 

fN C"' 

rT i 

A«,= (K-vJUUay, 

•XU. 


(B-38) 


(B-39) 


(B-40) 


(B-41) 


A(><i+i-k+t+')! t 
• CM>'. t'. Ci«,+ Jy)t-\+lt-i-»M 

«nd O^j Is definsd by BqusCion (fr>17). A synbol such ss dsnotss 

chs Issssr of the two qusntitlss inside the brecket. 




i+j+U+v -u-ir)>. 
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o)tf] 


(B-44) 


‘T- (B-45) 

D * \ (B- 46 ) 


For tha cace of t%io aquivalont eloctront and are Identical 


we have the following aliq)le expreaalone, 

where T and D are now given by 

O' s aUa-vK-vjp 


(B-47) 


(B- 48 ) 




(&- 49 ) 


and 4 (t,D) la given by Equation (B-44 ). 

A nuayber of formulae using the above equations are given in 
Table B-3. Machine computations for «nd have also been 

completed for all two electron combinations cheek marked in Table B-4« 
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TABLE B-3 


SLATER INTEGRALS 

(Th« formulae listed here are for F* (op; *nd G* (op; X^l), defined 

by Equation (b- 50), where subscript op are omitted for 

brevity.) 

Ci 4 '% 5 ^(hY% 

r6s,!^P)= \[\- 04xf(3)'4\)] 


FY3s,> 4; = i [ 1 - (\*X)%n 138 X+44X+'oy+| 

-ii[i - 04To\to)[%‘+5V4(&X+a)j 

F*(,p,34> - O+TO*’ (?5off 

^\0 8X417X+3;] 
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TABLE B-3 (continutd) 


F (3p,?JJ= -^^1 _ (3e6}'%|(,S'y\\\oX 

-floX-t\)] 

r 

F*C^p,i<^) = P “0+XJ (875 "X'h y‘47]l!V‘1\7X+4kiX*^ 


+-j?oX^iioV^aiy+a)| 


G-' 6 s,is)*a )t 0+xJ^(s 3f’‘- IS Xh. ( 3 ;) 


G'’0s,3s)= J x\i+):J^(s-x'^-5.X%t«X-iz.X4a-i) 
G'6 v|>)= |h x\i+x)^ 

G 6 s>fc):r y X^(l+Xy^(5-fcX- 8HX4 1-0 
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TABLE B-3 (continued) 


^4■3oX^>?9Sy->^83X-^t^ 

3 i -5 

G(is,3c|J ^ is X 

a X^(HX)"(ii‘+X-^.^X+9; 

G'(3s,}|)= i. x;'tw7X**-'J5-98 
QVir,3s)* A y^L\ti) 

G^(»|p.3p> ^-y^O^X) (ux^'M-s-y+jffJ) 
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TABLE B-3 (contlnuad) 


\74 C'+>^3 


\ ^ V 

G\Jlf>,^A> 1^ X^(t+r) (2^o5■)'-^l«?^-^UtJ) 
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The set of simultaneous equations (b-20) contain partial derl- 
k k 

vatlves of F (a,p) and G (a,3) with respect as well 

k k 

though F (a,^) and G are symmetric In the two electrons cx and 

P, their partlals are not symmetric. From the preceedlng expressions 
k k 

for F and G the following general properties of these two Integrals 
can be derived very simply, 

where R stands for either Integral F or G and 

The following expression for the partial derivatives are also derived 
easily, 

= R*^6^P>T»<f9 (1-55) 

T)1U^ 
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It Is therefore clear that both partial derivatives of F and 
G Integrals can be written as functions of a single parameter 

The case for t%fo equivalent electrons Is a particularly simple 
one; the partial of r^(a,a) la merely a nusdier, i^lch la. In feet, the 
value of this Integral fo*'^ » t*'® case of hydrogen atom (with 

atomic number unity) . 

=< ■ THE SCREENING PARAMETERS 


We now rewrite Equation (6-20) explicitly In terms of the 
partial derivatives of the Slater Integrals 






Is merely a number, Independent of any of the 


56) 

(3-57) 


(B-58) 


(B-59) 


f^C'i" SL) and g^C'i" SL) denote the coefficients In which the SL depen¬ 
dent contributions to the electrostatic energy are explicitly included. 

The syadx>l "i" denotes the orbital quantum nusters of all the incompletely 
filled stalls Including the parentages, etc. 
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(B-60) 


I 


He note thet. 


%(ft - 




and for the caee Z -» oo 


y . 'Y°* 


(B-6I) 


The nethod for coiq>uting the ecreenlng parametere ie aa followe. 
"a 

He flrat calculate a aet of n'a ulth X . = . He call theae aayoptotic 

Op 

n'a or . Hlth theae 's, and cone auitably chooaen large value of 
Z, we compute a new aet of and uae theae to recalculate x^'a by 
Iteration. He then decrease Z by aucceasive atepa and keep repeating the 
same Iterative procedure. Thus we work along an lao-electronlc sequence 
and obtain the screening paraiaeters for any desired number of atoms and 
Ions belonging to this Iso•electronic sequence. Since the procedure Is 
very time consuming, a program for computations on a digital computor 
has been developed* 


The asymptotic screening constants for two electron combina¬ 
tions, and using the average energy of a configuration only, can however 
be computed quite easily. They are defined as follows: 

(B-62) 

for a 4 
and 
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<B-63) 



The asymptotic screening constants for the average energy of a 
configuration can then be constructed from these two electron screening 
constants as follows: 




(B-64) 


b.8 the calculation of the transition integrals 


The only case of interest for an electric dipole transition is 
when the azimuthal quantum number of the active electron changes by one 
unit. Therefore, the calculation of the following integral is of suf¬ 
ficient generality: 


= 9 , 9 . 


(B-65) 


Jo 




(B-66) 




A similar integral with been evaluated by Gordoii^^^^\ 

The necessary formulas for the evaluation of Integrals involving 
confluent hypergeometric functions have been given, among other, by 
Landau and Lifshits^^^^^ in their mathematical aRpeadix. 
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The p index cen be reduced by using the foroule 


67) 


wherees the s Index cen be reduced by using another foroule . 

^ s(V -1 + S+aw»)Ty(- rt,t-Jrti sTjr (- V>rr«u4') 

Successive applications of Equations <6*67) and (6-68) to (6-66) yield 
•p\f2 / if- 4 Cp ^ 

(-VV^ ,-Vli n, J (f+Wt)J 

• Ji? (- ‘''•V Y'^ 

-tFri 

jf'f can be evaluated in closed form by another forsula 
ao -«la V' y -rfl— 

3;(«i.,*v),c-« ‘ a'RrXV^+V.) (Y»->0 (..„) 

(‘'..•'or'i- 4Y, V^V^-W*) 

where 2^^ « hypergeometrlc function given, In general, 

c/«k...,in . u a.Ca40t(V+0 

It clif^Tc^^O~r 

^ 0^ Cq’»»Xa4a) b(i>4iXb^a) 3^. 

c(c+i)(c+k; 3'. ' 
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The series termlnetes after a finite nusiber of terms irtien e end/or b are 
negetive Integers. On substitution of (B-70) in (B-69) end sinplifice- 
tion we get ^ p 

fM'.i <»•->'> 

where 

CT^”ri-> 


(B-72) 


(B-73) 


The following relation nay be used to reduce the nuster of hypergeometric 
functions from three to two in <B-72)9if desired. 


(e - aa, - cw- -» C*-V, 


(B-7*) 


wn«r« 

Combining (B-6S) with (B-72) we get 


&.-V,) 




-^IfO 






(B-75) 


(-Wf-o 
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which for th« spocUl rodueoi to tho following 

•quatlon first dorivod by Gordon^ 

tk ^ ^ 


(«r*0 




tt-76) 


Boocial caaoB can ba obtainad aaailv 


Vl-S 


Tha following apacial casas can ba obtainad aaaily 

li(*'p.«vy;V'*.)=3^ Iff ®- 

T . Iv'fJ'iV.Vt) -■JiCw-'W'vr ["Vr (»»■**•) 


■77) 








M-H 




T.(af7»'‘ >%y^) -C-^r'-^v 


(B-78) 


CB-79) 


(B-tO) 


- 3 

CY.->.r 

Orb)”-'** 
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B.9 ALTERNATE EXPRESSION FOR TRANSITION INTE(HtALS 


For th« sake of gonarality we shall evaluate the following 
Integral, 50 


“ .r** 


(B-81) 


For Coulood) potential wave funetlonst we shall later put q equal to aero, 
p equals unity for the electric dipole case, and It equals 2, 3, 4, etc. 
for higher sultlpole cases. Substituting from Equation (3-12) we get. 


(3-82) 


This Integration can be perforaed using certain standard results In the 
theory of confluent hypergeosatrlc functions, (see for exaegle Slater^*’^*\ 
pegs 54). isrfo{nia| the integration we obtain,. 


P«gs ^/. wcforniae kr* wm wulsw,. . J 

l5 , 9 , 9 , i ; . - 

• • M-***i>l W-<. 

• (B-83)* 


4A slsdlar result has also been obtained by Prof. D. H. Msnsel (wgniblished). 
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where the expression in the squere brecket following the susMstlon sign 
cen be readily seen to be equal to hypergeoaetrie func< 

tlon defined by Equation (^-71). The computation of the expression on 
the right hand side of Che above equation Is quite straight forward In 
principle, but becoass very laborious as n^ and n^ become large. 


He shall now proceed on to the facial case we are Interested 


In, namely the case where p - 1 and q > 0. Fbr this case, 

tF, (-Mi,{|+-Pi+4-+i i 

. J Ow?)l MtAO! 

i+j)'. i'. ~ 

= li 4 /n«-Y 

ii'.y,Vv) = ^, «»v i^'*^'***’’ 

>1"'^ 


(B-84) 


(B-85) 


|r £ £ A ii >*■* 

V >i 
' (vv3^ 


<B-«6) 
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to say that tha turn axprasslons can also ba obtainad by intagra- 
tlon of Equation (b~ 81) with tha radial wava functions givan by Equation 
(B-12) coad)inad with EquationCB-lA). Sons furthar siaplifieations of tha 
above expression can be brought about as follows. The conson demndnator 
of all tha tarM undar tha two suasuitioiis is ^ Tdking 

tha cosawn danoodnator and perforsding binonial axpansion of the resulting 
t.TBt lUMly <|5 +y 2 )'*^ obtain 

yHf+vf’ JiT , 






7^ 


W JK IMlHWi 


*' ■i- i-rt ip'. wi >' 




(B-87) 


laplacing p by another index u, such that u > p 4 J, wa gat, 

^“*3' 






Uso 




(B-88) 
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.r. 4.ci?i 


* Cv'f* < -W)l C 1)I 


u takes on all integral values free 0 to (n| + nj). f Is a funetlon of 
the quantum numbers of the two electrons under consideration, but it is 
independent of the screened nuclear charges, andy^ 2 * 
n^ are so large that desk calculations are too tedius, the calculation 
of on electronic computors is a trivial matter, giving us analyti¬ 
cal expressions of in terms of as follows, 

. »l (B-M) 




W^4«?’-a w 


rT'r 

^ 11- 


electrons, t is not. 


Table 2-5 gives expressions for Z^ for a seleeted number of 


transitions. 


tahj: B-5 dxpoli iniigrals 



b¥3tf±.. 

- 

I (3p,U) - 


-h'h) 

I,(4f» - 


1 

• 

1 

I,(5f.U) - 



X(6f.l-) * 

. 

W 

1 

1 
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TAMJ B-5 (eont.) 


1 

I 

I 

I 

I 


I,ar.u) - 

Sfx ^fx, 

I,(3f,2.) - [-3y, ^ ^5 


i,(tf2^ = ^ g?ty 


jjiei^'-*' 






121 



TAM.E B-5 (cont.) 




I 

I -- 

I 

T = 

- 
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TABLEB-5 (eont.) 


+2^4(6^‘J^ 

t224(/1 

-'00«(VK"ojJ|Ni‘>ji^l)(^.’a‘J 

.«20(3^Y^ -78«0/(fj 
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TABLE B-5 (eont.) 
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TAK.SB-5 (cent.) 






- 1 + 764 ^, yj 
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- 


I 

I 

1 


TAK.B B -5 (cont.) 


- 14 (^ ’ y}?f 

-Z«.4(3(g‘jNj.y)(S.’3.) 

•t59>5b^?^f| 
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TAMtS b> 5 (eont.) 


.Zt{Z5^ t 
.32(3SjJ 1350^^ 
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l(3i,Z^) • 


liWf) - 

l(7J,Zf) = 

m^) - 




TAS.E B-S (eont.) 
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Wf) - 
I(6<l,4f) - 


I/J r ^ ^ 

zsscj) y,X 

(rn.r ^ ' 

w" V* ''ft 

(Y,+YO* 

i»s6Cs*?Yi'^ 

CV.^tuV*- 


izeo(^h*j» 










TABLE B-5 (eont.) 
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TAM.BB-5 (eont.) 


I|«.y)- MVv) 

iw.ij)- ^q i g:gf (^j.-» 
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TAILS B*5 (eont.) 


liUM) * 

' 521 ^^ 

K^'-w) = [-S’ 

Ifh lyiL 

■ln,-’*>.i(.*''».i.-*^*l 
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K5{^j) 





B*5 (cent.) 


-‘f 


BIO Results of Conmutations ' 


(a) The Normalization Factor, 


is defined by Equation (B-13) 



1 _ 

(Yi+1): 


(B-13) 


The following recursion relations are useful in extending a 

table of 

ll/2 

'*n,i-l 


Q ^ Q ^ ^ (B-91) 




r " 1^/2 

sjI a+iii 

^ni |_ n n-i J '^n-l,f 


A few special eases of are given below: 

,-l/2 






i-l 

(B-93) 


(B-94) 


‘in.n.l • (2n)! 


1-1/2 


(B-96). 


A simple procedure In calculating Is to start with which 
• 1/2 

is always 2 , and successively calculate Q^ 2 > using Equation (B-91.) 

until one arrives at <1^ Equation (B-96) can be used as a final check 

* The work perfoneed for this section was sponsored by both this contract 
and Contract Ro. AF 33(600)-42617. 
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The values of for n up to 9 and all allowed values of I were 
computed in the above isanner, and subsequently machine confutations 
were also performed. Table B-6 lists the 

(b) The Coefficients d|^^ 

The d|^^*s are the coefficients occurring in the hypergeosmtric 
function (-(n-i-1). 2i-f2:x^ . They are given by Equation (B-15). 


“■n/ ■ <■» (.Sllllii (jHltlii iV 

(B-15) 

For 1 > n-f-1, d^ vaalshat. 


The following recursion relations are useful 
the table of d^: 

in extending 

.i fn-i-i) .1-1 

‘*ni • (21+1+1)1 **01 

(B-97) 

A . _(n»l-i)^(21+l) i 

“nl (n-lH21+i)(21+l+i) “n.l-l 

(B-98) 


(B-99) 

‘*nl * (2i+l{^+l+l) ‘*n-l.l-l 

(B-lOO) 

A few special cases of d^ are of interest 



(B-101) 


13S 


. 2041241 ^ 


.93S414»^ 


.3333331+00 8d .3622841+00 

.372678B-01 8f .639708B-01 

.4564358+00 8g .6155598-02 

.7905698-01 8h .3494708-03 

.4980128-02 81 .1185408-04 

.5773508+00 8k .2186208-06 

.1322888+00 9p .1054098+01 

.1259888-01 9d .4624818+00 

.5249508-03 9£ .9343538-01 

.6972168+00 9g .1046258-01 

.1972038+00 9h .7117658-03 

.2439758-01 91 .3060688-04 

.1515408-02 9k .8244688-06 

.4569118-04 9i .1249778-07 

.8164968+00 
.2738618+00 
.4123938-01 
.3290308-02 
.1465378-03 
.3386858-05 
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(B-102) 


n-X-1 
■ 2 (^+ 1 ) 


d 


1 

n,0 


(-1)^ 



(B-103) 


A table of for values of n up to 5 was prepared using the above 
recusion relations and subsequently machine computations were perforond 
for n up to 9. These are given in Table B-7 

(c) The Screening Parameters 

Equations (B-56) define a system of simultaneous non*linear 
equations in which the right hand side is a function of all the screening 
parameters of a given problem. A Fortran program for use 

on digital computers has been developed for solving these equations. 

Ue work along an isoelectronlc sequence. We use as our starting point 
the case when Z -* oo. The screening parameters for this case can be 
computed easily with given by Equation (B-61). As an intermediate 
step we perform iterative calculations for Z * where Z^ is the 

lowest msad)er of the isoelectronlc sequence, using n^'s for Z > oo case 
as our starting point. The iterations arc carried until the relative 
difference between the results of two successive iteration is less then 
a prescribed quantity e (for computations reported here € * 0.0001). 

The process is repeated successively for the cases Z ■ Z^ + 19, 18, 

....Z^. The screening paramsters, expressed as s(a)(* ^ n^) are given 
in TAbles B-8 to Bol9. 
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.7142KI-01 -.2976191-02 
,4166671-01 -.130M91-02 
.2727271-01 -.7575761-03 
.1923081-01 -.4578751-03 







TABLE B-7 (continued) 




S^SSSoooSSooSo^ 

(nd«(<4ro«r4nmomO>e«r>> 

«n»o«OMvOinirt<-)ioeover4<2Q 


00 ^ i-i lA '4' m 


??????????????? 

HHHOaMHHHHBdtOUHtOM 

r^«n«0<oO^^*n<^<^«»>r>*oo^»^r<.c»> 

«<n^mg'0«A>-ieo(o<0(nm«Ofn 

<er>io<Ae<eiA««co«ooi«o«m 

««A^«O2J«O«»»f«0e»>»OlAlA24«*> 

«e<n•-•ano^<oot(ncn«oo^s^e«(*> 

•-4(nr»iA(A'400Mf-iaOiHlAi-lf>4M 


ilsssiiails sill 










Ion 

TABLE B-8 

SCIBBMING PAXAMBTBBS 

Be I laoelaetronlc Sequence 

Configuration: l8^2s^ 
Term: ^S 

B(lB) 

b(28) 

Ll‘ 

3.2761 E-1 


1.7902 EfO 

Be I 

3.1861 8-1 


1.8599 B+0 

B II 

3.1369 E-1 


1.8917 EfO 

C III 

3.1184 B-1 


1.9025 EfO 

M IV 

3.1164 E-l 


1.9036 EfK) 

0 V 

3.1221 B-1 


1.9006 Ef-0 

P VI 

3.1312 E-1 


1.8961 EfO 

He VII 

3.1414 B-1 


1.8912 B+0 

Na VIII 

3.1517 E-1 


1.8863 EfO 

Mg IX 

3.1616 B-1 


1.8817 EfO 

A1 X 

3.1708 B-1 


1.8775 EfO 

81 XI 

3.1794 B-1 


1.8736 EfO 

P XII 

3.1873 B-1 


1.8700 EfO 

8 XIII 

3.1946 E-1 


1.8668 EfO 

Cl XIV 

3.2013 E-1 


1.8639 EfO 

A XV 

3.2074 E-1 


1.8612 EfO 

K XVI 

3.2130 B-1 


1.8587 BfO 

Ca XVII 

3.2182 B-1 


1.8565 BfO 

8e XVIII 

3.2230 E-1 


1.8544 BfO 

T1 XIX 

3.2274 E-1 


1.8525 BfO 

Z " 00 

3.3353 B-1 


1.8078 BfO 


140 




TA1W.lt B-9 

SCUBNIMG PARAMBTBIS 
BeT lAo«l«ctronic Sequence 
2 

Conflguretion: Is 2s2p 
Tens: 


Ion 8(Is) s(2s) s(2p) 


Li’ 

3.1832 

1-1 

1.5178 

BfO 

3.0030 

BfO 

Be I 

3.2321 

B-1 

2.1229 

BfO 

2.0640 

BfO 

B II 

3.2444 

B-1 

2.1965 

BfO 

1.9749 

BfO 

C III 

3.2590 

B-1 

2.2165 

BfO 

1.9390 

BfO 

N IV 

3.2766 

B-1 

2.2186 

BH) 

1.9180 

BfO 

0 V 

3.2949 

E-1 

2.2138 

BfO 

1.9037 

BfO 

P VI 

3.3126 

B-1 

2.2065 

BfO 

1.8932 

BfO 

lie VII 

3.3289 

B-1 

2.1987 

BfO 

1.8849 

BfO 

He VIII 

3.3437 

E-1 

2.1911 

BfO 

1.8783 

BfO 

Mt nc 

3.3570 

B-1 

2.1840 

BfO 

1.8728 

BfO 

A1 X 

3.3690 

B-1 

2.1775 

BfO 

1.8682 

BfO 

81 XI 

3.3797 

B-1 

2.1715 

BfO 

1.8642 

BfO 

P XII 

3.3893 

B-1 

2.1662 

BfO 

1.8608 

BfO 

8 XIII 

3.3980 

B-1 

2.1614 

BfO 

1.8578 

BfO 

Cl XIV 

3.4059 

B-1 

2.1570 

BfO 

1.8551 

BfO 

A XV 

3.4130 

B-1 

2.1530 

BfO 

1.8528 

BfO 

K XVI 

3.4195 

B-1 

2.1494 

BfO 

1.8507 

BfO 

Ga XVII 

3.4254 

B-1 

2.1441 

BfO 

1.8488 

BfO 

8c 8VIII 

3.4308 

B-1 

2.1430 

BfO 

1.8471 

BfO 

T1 XIX 

3.4338 

B-1 

2.1402 

BfO 

1.8436 

BfO 

8 

s 

N 

3.5M7 

B-l 

2.0773 

BfO 

1.8129 

BfO 
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TAH.E B-10 
8CREBMIIIG PASAMBTBBS 
B« I IsoaUetronlc S«qu«nc« 
2 

Configuration: It 2t 3p 
T*ini! 


Ion 

• (!•) 

• (2p) 

8(3p) 

Ll" 

3.1396 E-1 

1.9799 BfO 

2.5502 BfO 

B« Z 

3.1969 S-1 

1.8484 BH) 

2.7110 BfO 

B II 

3.2529 B>1 

1.7602 BfO 

2.7775 BfO 

C III 

3.2957 B-1 

1.7074 BfO 

2.7953 BfO 

N IV 

3.3280 B-1 

1.6774 BfO 

2.7864 BfO 

0 V 

3.3527 B-1 

1.6598 BfO 

2.7677 BfO 

P VI 

3.3720 B-1 

1.6490 BfO 

2.7471 BfO 

Me VII 

3.3876 B-1 

1.6419 BfO 

2.7276 BfO 

If« VIIZ 

3.4003 B-1 

1.6370 BfO 

2.7102 BfO 

Mg IX 

3.4110 B-1 

1.6334 BfO 

2.6950 BfO 

A1 X 

3.4200 B-1 

1.6307 BfO 

2.6818 BfO 

81 XI 

3.4276 B-1 

1.6285 BfO 

2.6703 BfO 

P XII 

3.4343 B-1 

1.6268 BfO 

2.6603 BfO 

8 XIZI 

3.4401 B-1 

1.6253 BfO 

2.6514 BfO 

Cl XIV 

3.4453 B-1 

1.6240 BfO 

2.6437 BfO 

A XV 

3.4498 B-1 

1.6230 BfO 

2.6368 BfO 

K XVI 

3.4539 B-1 

1.6221 BfO 

2.6306 BfO 

C« XVII 

3.4575 B-1 

1.6212 BfO 

2.6251 BfO 

8c XVIII 

3.4609 B-1 

1.6205 BfO 

2.6201 BfO 

Ti XIX 

3.4639 B-1 

1.6199 BfO 

2.6156 BfO 

X ■ 00 

3.5261 B-1 

1.6072 BfO 

2.5266 BfO 
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MKJ B-ll 

t(»imii6 rAIAMETBS 
0 I iMMlMtronie S«qu«iic« 
2 2 4 

ConfifurAelon: If 2fl 2p 
T«ni: ^8 


Ion 

• (!•) 

«(2f) 

• (2p) 

H’ 

3.3288 B-1 

2.9181 BfO 

3.6713 BfO 

0 I 

3.4078 1-1 

3.0238 IfO 

3.5794 BfO 

F II 

3.4825 B-1 

3.0779 BH) 

3.5248 BfO 

Hf III 

3.5524 B-1 

3.1080 BH) 

3.4873 BfO 

Hf IV 

3.6170 B-1 

3.1253 BfO 

3.4595 BfO 

Hg V 

3.6763 B-1 

3.1354 BfO 

3.4378 BfO 

A1 VI 

3.7303 B-1 

3.1411 BfO 

3.4204 BfO 

Si VII 

3.7794 B-1 

3.1442 BfO 

3.4061 BfO 

P VIII 

3.8241 B-1 

3.1455 BfO 

3.3940 BfO 

8 DC 

3.8648 B-1 

3.1458 BfO 

3.3837 BfO 

31 X 

3.9019 B-1 

3.1453 BfO 

3.3748 BfO 

A XI 

3.9358 B-1 

3.1445 BfO 

3.3671 BfO 

K XII 

3.9668 B-1 

3.1433 BfO 

3.3602 BfO 

Ca XIII 

3.9953 B-1 

3.1420 BfO 

3.3541 BfO 

8e XIV 

4.0216 B-1 

3.1406 BfO 

3.3487 BfO 

T1 XV 

4.0458 B-1 

3.1392 BfO 

3.3438 BfO 

V XVI 

4.0682 B-1 

3.1377 BfO 

3.3394 BfO 

Cr XVII 

4.0890 B-1 

3.1362 BfO 

3.3353 BfO 

Mn XVIII 

4.1083 B-1 

3.1348 BfO 

3.3316 BfO 

P« IDC 

4.1263 B-l 

3.1334 IfO 

3.3283 BfO 

Z • 00 

4.6176 B-1 

3.0813 BfO 

3.2504 BfO 
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The hydrogenlc wave functions using the screening parameters 
are compared with the self‘consistent field wave functions for neutral 
oxygen (Hartree, Hartree and Swlrles^®”^®^) and negative Ion of lithium 
(Roothaan, Sachs, and WeissIn Figures B-1 to B-4. The hydro¬ 
genlc wave functions are given below: 

0 I (ls^2e^2p^ ^S) Is electron: 42.394 r exp (-7.6592 r) 

0 I (ls^2aV ^S) 2s electron: 7.8491 r exp (-2.4881r) (l-2.4881r) 

0 I (ls^2sV ^S) 2p electron: 8.3864 r^ exp (-2.2103 r) 

LI (l8^2s^) 2s electron: 0.94079 r exp (-0.66524r) (l-0.66524r) 

The agreesmnt between the two sets of wave functions Is 
resmrkable for 0 1 and la good for LI. 

Total Ionization energies based on hydrogenlc wave functions 
are compared with experlsKncal values In Table B-20. The agreement Is 
excellent; for Be I It Is as good as that for the multlconflguratlon 
Hartree-Fock wave functions of Klbartas, Kavetskls and lutslsf^ 
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Figinrc 1-4. Li" Wave Function of 2s Electron 


TABLE B>20 


TOTAL IONIZATION ENERGIES (Atomic Unit) 



Ifydrogenlc 

Wave Function 

H-F 

Wave Function 

Experimental 

0 I (4) 

- 74.62 


- 74.95 

P II (^S) 

- 98.61 


> 98.92 

Be I (^8) 

- 14.70 

- 14.64* 

- 14.67 

B II <4) 

- 24.38 


- 24.35 

C III (^8) 

- 36.55 


• 36.55 

N IV (^8) 

- 51.23 


• 51.25 

0 V (^8) 

• 68.40 


- 68.45 


*BaMd upon multiconf^urftlon Hartrcc-Fock calculations of Klbartas, 
Kavatskla and lutalaJ®"^*' 





(d) The DipoI« Integrals 


The electric dipole integral is given Equation (B*90) 
in terme of a coefficientY'^Cnj^^^* defined by Equation (B»89). 

Coaputations on an IBM 7090 have been perfornsd for the coefficients 
T” for all two-electron coabinatione of practical interest* ai^ the 
results are given in Table B-21 to B-26. In order to reduce huoan 
errors in the reproduction of a large collection of nunbers, the print¬ 
out froa the digital coaputor has been photocopied. P8I 7* for exaaple* 
aeanst^^. Nuadiers such as 0.23999999 or 0.24000011 should be rounded 
off to 0.24. 

2 

A few saaple values of the transition integral o coaputed 
with the screened hydrogenic wave function are given in Table B-27* 
and compared with calculations based on Coulosd) approxiaation^^"^^^ 
and H-F wave functions. The agreeasnt is satisfeetory. 


U-t 



159 






























160 

























XUU >>11^ (kjS.AjP) (centloMd) 



161 


rSl $»C.I»2iliS976ki 































B2P> (eoBtlaiMd) 




162 











































163 

















164 



















165 




































lABLB B-23 ^ (n^P.n2d) (eontimiMi} 



167 














16S 



































(4>7) (coacinied) 



169 













170 













171 






























TABLE B-27 


TRANSITION INTEGRAL cr^ 


Atom 

2 

Transition Integral o' 

Ifydrogenlc 

Coulomb 
Approx. 

Hartree-Foek 

Be I 

2.27 

2.03 

1.86 

B II 

0.977 

1.010 


C III 

0.544 

0.573 


N IV 

0.348 

0.368 


P XII 

0.0522 

0.0540 
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APPENDIX C 
D. Layser 


NEW GENERALIZATION OF THE SCREENING THEORY 
AND THE HARTREE-FOCK THEORY 

C .l. Introduction 


The distinguishing feature of the Z-dependent theory of 
atomic spectra (Layzer/^ henceforth referred to as A) 

is the fact that it allows for certain effects of configu¬ 
ration mixing in the first approximation. This ensures that 
the predicted term energies have the correct linear depen¬ 
dence on Z. Experiment shows that the structure of term 
spectra, in both simple and complex systems,, depends strong¬ 
ly on the linear Z-dependence of the energies. A knowledge 
of the energy coefficients and in the expansions 

W = + WjZ + Wq + 0(z‘^) (C-l) 

enables one to describe quantitatively the most conspicuous 
features of term spectra in isoelectronic sequences (see 
especially Figures I and 2 in reference C-l). 

In order to calculate the coefficients (W^ is given 
by the well-known hydrogenic formula) one needs to know 
nothing about the radial dependence of the atomic charge 
distribution beyond what is provided by the theory of the 
hydrogen atom. That is, all the necessary calculations 
utilize hydrogenic wavefunctions. One can, of course, use 

less simple wavefunctions in these calculations, as was 

IC-i) 

done in the first calculations of this type (Layzei ' 
but such refinements only improve the accuracy of higher 
terms in the expansion. 
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For some purposes, however, one does require a reason¬ 
ably accurate knowledge of the radial charge distribution 
and hence of the radial wavefunctions. One needs this in¬ 
formation, for example, in order to calculate electric- 
dipole transition probabilities. 

The various procedures in common use for calculating 
accurate radial wavefunctions are all based on the varia¬ 
tion principle. They differ through the restrictions they 
impose on the form of the variational functions. Hartree- 
Fock wavefunctions (i.e., exact solutions of the Hartree- 
Fock integro-differential equations) fall at one end of 
the gamut of variational wavefunctions, screened hydro- 
genic wavefunctions at the other. The latter involve only 
a single adjustable parameter, while the former are re¬ 
stricted beforehand only by the usual conditions of inte- 
grability and smoothness. Other types of variational wave- 
functions seek to achieve some of the accuracy of Hartree- 
Fock wavefunctions without sacrificing all the simplicity 
of screened hydrogenic wavefunctions. Whicli type of varia¬ 
tional wavefunction one chooses to employ will be governed 
by the purpose in view. The Hartree-Fock description is 
certainly the most accurate; from a qualitative standpoint, 
however, the screening description offers unrivalled ad¬ 
vantages of simplicity and flexibility. 

Any meaningful comparison between different variational 
wavefunctions presupposes a common approximation scheme. 
Thus we may compare wavefunctions calculated with neglect 
of configuration mixing, or with certain specified kinds 
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of configuration mixing taken into account. (In practice« 
the calculations that allow most fully for configuration 
mixing are those based on the simplest wavefunctions, so 
that strict comparisons are difficult.) 

According to the Z>dependent theory, the most is^rtant 
kind of configuration mixing is that between configurations 
belonging to the same complex (i.e., the same set of princi- 
pal quantum numbers and the same parity). On the other hand, 
the theory of radial wavefunctions has usually been devel¬ 
oped within the framework of the central-field approximation, 
though attempts have occasionally been made to allow for 
effects of configuration mixing. A notable example is the 
work of Hartree, Hartree, and Swirles^^'^\ which dealt 
with the mixing of configurations ls^2s^2p” with configura¬ 
tions l8^2p”^^. However, such attempts have lacked an ade¬ 
quate theoretical framework, comparable in simplicity and 
generality with that provided by the central-field approxi¬ 
mation. The present paper seeks to supply such a framework. 

C .2. The Theoretical Framework 

Let Y denote an antisymmetric N-electron wavefunction 
pertaining to a definite state in a definite electronic con¬ 
figuration. The quantity 

W = <YHY> , (C-2) 

where H is the N-electron Hamiltonian and the brackets <C y 
indicate an integration over all 3N coordinates and a sum¬ 
mation over all the spin variables, is the variational ener- 
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gy. It is a functional of the radial wavefunctions. If 
these are specified beforehand, except for the values of 
certain parameters, then W reduces to a function of these 
parameters. W exceeds the true energy of the state under 
consideration provided that ^ is normalized and is orthogo¬ 
nal to the true eigenfunctions of all states lower than the 
one under consideration. The condition 6W » 0 leads to a 
set of simultaneous integro-differential equations for the 
radial wavefunctions, or to a set of simultaneous algebraic 
equations for the adjustable parameters. 

One can easily generalize the theory just sketched so 
as to include the effects of mixing among a finite number 
of configurations. Let denote variational 

wavefunctions for states with fixed values of the quantum 
numbers SLp in different configurations. Suppose, for the 
sake of definiteness, that corresponds to the deepest of 
the states under consideration. Then the corresponding 
variational energy is the lowest eigenvalue of the matrix 
II II . As before, the condition 6W^ » 0 yields a 

set of simultaneous equations for the unknown functions or 
parameters. (Note that all the radial wavefunctions or 
parameters that figure in the complete set may 

occur in a single equation.) 

An accurate solution of the secular equation for 
would lead to a set of variational equations of prohibitive 
complexity. The procedure usually adopted consists in 
retaining only first- and second-order contributions to 
in the sense of ordinary perturbation theory. Thus is 
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given by 




r 




(C-3) 


where denotes the mean energy of the ith configuration. 
Even with thia sifl^lificatlon, the variational equations with 
configuration mixing taken into account present considerable 
more difficulty than the usual equations. And frequently the 
effects of mixing are so large that second-order perturbation 
theory does not describe them adequately. 


Another important difficulty inherent in the usual methods 
of dealing with configuration mixing concerns the choice of 
interacting configurations. It is of course true that the 
more interactions one takes into account« the more accurate 
will be the theoretical description. It has been amply dem¬ 
onstrated, however, that in order to obtain predictions of 
term spacing in neutral and slightly ionized systems that 
are accurate to within a few percent, one needs to allow for 
a very large number of interactions. Hence a calculation 
that allows for mixing between a small number of configura¬ 
tions would be significant only if the choice of the configu¬ 
rations considered could be justified on grounds other than 
ma thema tica1 convenience. 


Hie problem of configuration mixing takes on a differ¬ 
ent complexion if we adopt the Z-dependent approach. Con¬ 
sider the expansion 
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(C-4) 


*= Evi' 

i 

where the refer to distinct configurations and the sum 
runs over all ccmfigurations connected with configuration 1 
(say) by nonvanishing energy matrix elements, (we assume 
that the state corresponding to has been nominally assigned 
to configuration 1). 

The numerical values of the coefficients a^ depend on 
the choice of radial wavefunctions. For example, if the 
radial wavefunctions in are solutions of the Hartree- 
Fock equations, then the coefficients a^ pertaining to con¬ 
figurations with the same set of azimuthal quantum numbers 
as configuration 1 will all be very small, because the vari¬ 
ation principle ensures that most of the effects of mixing 
between configurations characterized by the same set of azi¬ 
muthal quantum numbers have already been incorporated in the 

wavefunctions. Now, every coefficient a. is either 0(1) or 
-1 ^ 

0(Z ). The coefficients of the first type are associated 

with configurations that have the same set of principal 
quantum numbers as configuration 1; those of the second type 
with configurations that have a different set of principal 
quantum numbers. 

The coefficients of the first type can be written in 
the form 

ai = a° + 0(z'^) , (C-5) 

0 

where a^ is independent of Z. It was shown in A that the 
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constants do not depend on the form of the radial wave - 

functions . The coefficients and parts of coefficients that 
do depend on the form of the radial wavefunctions are 0(Z~^). 

We may conveniently split off the terms on the right side 
of ( C-4) that have constant coefficientst 

i i 

where 

a^ = a^ - a® » 0(z“^) . (C-7) 

The nonvanishing coefficients are finite in number. They 

can be calculated with arbitrary precision, being essentially 
by-products of the calculation that yields the energy coeffi¬ 
cients (see A). 

We propose to replace the usual approximation, which con¬ 
sists in neglecting configuration interaction entirely, by an 
approximation that represents the wavefunction for a given 
state by the finite sum 


f - Y. “I'^i ’ 

where the coefficients may be regarded as Known numbers. (For 
spectroscopic terms belmtging to completes of the type (1^2^) 
and (1*2^)° they are given in Table C-i. 


The approximation just described becomes increasingly 
accurate as z increases along an IsoelectrMtic sequence. 


183 



184 



Even for moderately ionized and un-ionized systems, however, 
it should prove superior to the central-field approximation, 
since experiment has shown that the linear Z-dependence of 
the term energies usually dominates the structure of the term 
spectrum. 

C. 3. The Generalized Screening Parameters 

Notation. The letters r,t,t' will denote configu¬ 

rations; the symbol P will distinguish different states in the 
same configuration; a will denote an orbital nl; and a will 
denote a complex. Thus the label a(r)r distinguishes the 
state r nominally assigned to the configuration r of the com¬ 
plex a. In the present notation, equation ( C-4) takes the 
form 


f(arr) = ^<arrltr>'i'(tr) . (C-9) 

t 

The Hamiltonian H for a many-electron system is given, 
in the nonrelativistic approximation, by 

H = IIq + V , (C-10) 

where 


Hq = ^ (-Pi/2m - ZeVr^l , 


(C-li) 


V 


i<j 


(C-IZ) 
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The variational energy in the present approximation is 
defined as 

W = <arf|n|arr> . (C-13) 

We now evaluate W using screened hydrogenic wavefunctions 
for the undetermined radial wavefunctions that figure in 
this formula. 

For the one-electron part of W we obtain 
<arrlHQlarr> = <arri tr> < tPl H^l t • r>< t * Pj arP) 

t, t * 

= Y l<arPltP>l2 <tPlHpltP> 
t 

t a (C-14) 

where denotes the number of orbitals of type a in the 
configuration t, and s^^ denotes the screening parameter 
that corresponds to the orbital a in configuration t. A 
more complete notation for this screening parameter is 
s^(at|rP). The suffix Z indicates that the screening para¬ 
meter depends in general on Z; s^ differs from by terms 
0(Z . The labels rl’ indicate that the value of the screen¬ 

ing parameter depends on the state under consideration. 


For the two-electron contribution to W we have 


I 


<arr|vlarr> = ^<arr|tr> <trlvlf r> <t'r|arr> . 

t,t' (C-15) 

Hence 

<»rrlHl»rO - ^|<atr|tr>l^ * 

t a 


The variational equations are 




(C-i6) 


3r-- - sr; «triv|to + ,rr<‘nvifi» - o. 

”o tyt (c- 17 ) 

where 

V ® - V - • 

(C-U) 

The diagonal matrix element ^tr|v|tr^ may ^ written in 
the form 


<tr|v|tr> - , ,c.i„ 

Where the contribution represents the average value of 

<tr|v|tr>. Explicitly, 


V 


t 


^<taSL|vltaSL>(2S + 1)(2L ^ 1) 
aSL 

(2S + 1)(2L ♦ 1) 

aSL 


(C-20) 
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The term-dependent contribution vanishes when the con¬ 
figuration t consists of complete groups only, or of complete 
groups plus or minus a single electron. The contribution to 
v^p from interactions between electrons in a group l” satis¬ 
fies the symmetry condition 

v(l"r) = . (C-21) 


In a first approximation one may wish to neglect the term- 
dependence of the screening parameters. We accordingly define 


3v 

a (at) = ^ X . (C-22) 

ot 

If the configuration t contains no more than one inccxnplete 
group, can be put in the form 


where 


''t ■ ' 

a^ 


(a = p) 



(a ji p) 


(C-23) 


(C-24) 


The two-electron interaction (ap) is defined in A; in a first 
approximation (aP) % F^(ap). The coefficient denotes the 

number of ways in which two electrons labelled by the indicated 
quantum numbers can be selcted from the configuration t. 


With the help of these definitions we can write s(at) in 


the form 


s(at) 


(C-25) 


” I % ■ ‘c.p>’<“lP> ' 

p 

where 

a (alp) = (1 - 2 (ap) • (C-26) 

The two-electron screening parameter s(a|p) provides an unam¬ 
biguous, qunatitative measure of the screening of an a elec¬ 
tron by a p electron. 

The screening parameters s (a|p) are independent of the 
configuration t, and can therefore be evaluated once and for 
all. 


The exact screening parameters are given by 
2 


s(at|rr) 




^5!‘^tf;rr''ttT 

tyt 


) , 

(C-27) 


where 


Vtt r " <tr|v|tT> . 


(C-28) 


In this equation one is not in general justified in neglecting 
the contribution from configuration mixing. In some instances the 
second term in the brackets ( ) may be considerably larger than 
the first term. This is always the case, for example, when the 

configuration t is composed of complete groups but contains one 

2 2 

or more incomplete shells, e.g.. Is 2s . Then v^j, ■ 0, but ft 0. 
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Ttie aet of equations ( C>27) for all orbitals a in the 
mixture ( C-9) determine the corresponding screening para¬ 
meters s^iotlrP). In spite of their complicated appearance« 
these equations are comparatively easy to solve numerically. 

I., practice one uses an iterative procedure, first evaluating 
the right sides for Z a This gives the screening constants 
s^(at|rr). In the next ai^roximation one evaluates the ri^t 
sides for finite values of Z, replacing the parameters s by 
the constants s , and so on. 

M 

One can take advantage of the dominant contribution of 
s(at) to s(at|rr) to shorten the iterative evaluation of Z- 
dependent screening parameters. One begins by solving the 
system of equations (C-22) for s^iat). One then evaluates 
the term-dependent correction to s^CatlrP) using s^Cat) in 
place of S 2 (at|rP) in equation ( C>27). with this method 
one needs to go through the iterative procedure only once 
for a given configuration. 


C.4t 9fpgyatj|.|;ed Hartygg-fock 


The Hartree-Fock equations may be written in the form 


( - i- 


1(1 ♦ 1) _ ^nl 
2r^ 




(C-29) 


The quantities and are functionals of all the radial 
wavefunctions appropriate to the configuraticm being consid¬ 
ered. The exchange term has its origin in the antisymmetry 
of many-electron wavefunctions. 


If we take as our starting point the a^roximation ( c-9) 


IfO 


instead of the central-field approximation^ we arrive at an 
equation that is formally identical with ( C-29) This equa¬ 
tion differs, in fact, from the ordinary Hartree-Fock equation 

only through the exchange term X ,, which contains additional 

n i 

contributions arising from configuration mixing. These con¬ 
tributions are of precisely the same type as the ordinary ex¬ 
change contributions. Thus the increase in accuracy that re¬ 
sults from allowing for the dominant effects of configuration 
mixing has not been bought at the price of added complexity. 

There exists a very close analogy between the Hartree- 
Fock theory and the hydrogenic screening theory, which it is 
convenient to exploit because the screening theory is so much 
simpler than the Hartree-Fock theory. The generalized Hartree- 
Fock equations are analogous to the screening equations {C-Z7). 
Just as in the screening theory one is not justified in neglect 
ing the mixing contributions of beside the term-dependent 

contributions so in the generalized Hartree-Fock theory 

one is not justified in neglecting the mixing contributions to 

X , , 

nl 

Calculations based on the generalized Hartree-Fock theory 
are planned for the near future. 
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APPENDIX D 

P, 

APPLICABILITY OF THE SCREENING THEORY TO HEAVY ATOMS 
D.I. Introduction 


In the following^ the term "screening theory" will be 
used to designate the Z-dependent approximation scheme descri¬ 
bed by Layzer and, in a generalized form that allows 

fully for relativistic effects, by Layzer and Bahcall 
The term "Slater theory" will be used to designate the con¬ 
ventional theory of atomic spectra, whose basic approximations 
were introduced by Slater in his classic paper of 1929p’^^later 
recent two-volume monograph, "Quantum Theory of Atomic Structure", 
describes the modern form of the theory. 

The screening theory was devised to account for certain 
systematic discrepancies between the Slater theory and experi¬ 
ment, pertaining to the centers of gravity of spectroscopic terms. 
These discrepancies are sometimes large enough to reverse the 
predicted order of term energies in a configuration. 

Now, any reasonably general theory of atomic structure mui-t 
incorporate some fairly drastic simplifying assumptions The 
Slater theory rests on the assumption that the effects of con¬ 
figuration mixing can be neglected in a first approximation. 

The screening theory, on the other hand, rests on the assumption 
that the energy of an atomic system can be expressed as a rapidly 
convergent power series in (Z - s)~ , where Z is the nuclear 
charge and s is a screening parameter. This assumption leads 
naturally to a theory that is in most respects simpler than the 
Slater theory, but which allows for certain effects of configura¬ 
tion mixing in the first approximation. 
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As yet, the screening theory has been applied only to 
light atoms — specifically to atoms with ten or fewer elec¬ 
trons. In these applications it has proved very successful, 
accounting in a quantitative way for features of term spectra 
that the Slater theory describes only semi-quantitatively or 
qualitatively. The question now arises whether the screening 
theory will also provide a useful description of the spectra 
of heavy atoms. The answer hinges on how rapidly expansions 
of the term energies in inverse powers of Z - s converge. 

The fact that Z - s increases with N (the nun4>er of atomic 
electrons) argues for the applicability of the theory to 
heavy atoms. On the other hand, one would also expect the 
coefficients in the energy expansion to increase with N, 

Since it is far from obvious to what extent the second effect 
will offset the first, the question cannot be decided with 
confidence on general grounds; specific calculations are 
required. 

The iron-group elements provide ideal subjects for such 
calculations. They are sufficiently heavy to afford a defini¬ 
tive test. At the same time, they are sufficiently light to 
enable relativistic contributions to the energy to be treated 
as perturbations. Finally, the spectra of many atoms and 
ions in the iron group have been well analysed experimentally. 

In evaluating the merits of a theoretical description, 
one must compare it not only with the data it sets out to 
describe but also with other available theoretical descriptions. 
During the past few years, attempts have been made to patch up 
the Slater theory, as it applies to the spectra of elements in 
the iron group, by adding to the theoretical energy formulae 
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an empirical correction term of the form aL(L -t- 1) . That a 
term of this type could account for a fairly large part of the 
discrepancy between the Slater theory and experiment, in spectra 
of the iron group, was first pointed out by Layzer and, 

independently, by Racah.^’'” ' The Slater theory, amended in 
this way, has since been used extensively by Trees, Rohrlich, 

and Racah and his pupils in semi-empirical analyses of term 

(r-7) 

spectra in the iron group. Racah has suggested that a 

theoretical basis for the L(L + 1) correction is afforded by 
the work of Bacher and Goudsmit.^^'®^ 

In the following, we shall re-examine the theoretical 
basis of the L(L + 1) correction. We shall also present the 
results of new calculations, based on the screening theory, 
for complexes with a partly filled M shell. We shall find 
that the modified Slater theory cannot account quantitatively 
for the effects of configuration mixing described by the screen¬ 
ing theory. We shall also find that the latter are of the same 
order of magnitude as the observed discrepancies. This suggests 
thivt i:hc screening theory will provide an adequate description 
of these spectra. Unfortunately, sufficient experimental data 
are not yet available to teat the theoretical calculations 
presented here, the complexes here studied having been >hosen 
for their comparative mathematical simplicity rather than for 
comparison with experiment. However, work now in progress 
should lead before long to a definitive experimental test of 
the screening theory as applied to atmns and ions in the iron 
group. 
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Theoretical Results for Complexes with a Partly-Filled 
M Shell * 

We have considered the following complexes: , 

,,2.8_2,0 ,,2-8.16, ,,2.8-16,0 _ .. , 

(1 2 3 ), (1 2 3 ), (1 2 3 ) . The configurations belong¬ 

ing to these complexes are listed in Table u>i. These are the 
simplest complexes with partly-filled M shells, other than 
those consisting of closed shells plus or minus a single elec¬ 
tron. Complexes containing from three to seven 3d electrons 
present a somewhat more difficult theoretical problem. 

The evaluation of the energy coefficients and the 
unitary matrices that diagonalize the energy submatrices 

^/ \ CT proceeds along the lines laid down in the reference 
fD-l) 

cited earlier (Layzer ). is the coefficient of Z 

in an expansion of the nonrelativistic part of the energy in 
inverse powers of Z; '^(n)pSL ^ finite matrix whose elements 
connect states labelled by the indicated quantum numbers 
is an eigenvalue of the radial integrals that figure 

in the matrix elements being evaluated with straight hydro- 
genic wavefunctions.] 

The results of the calculation are present in Tables D-2 thru 

D-9. Tables D-2 thru D5give the coefficients Wj^ as well as the 

values that would have been obtained for these coefficients 

if configuration interaction had been neglected (i.e , the 

values of Wj^ according to the Slater theory) . Tobies D-6 r’ ru D-9 

give the unitary m.'trices that diagonalize the various energy 

subroatrices V, , . 

(n)pSL 

* The work summarized in this subsection was supported in part 
by a grant from the National Science Foundation. 
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TABLE D-1 


TBIM8 AND CONPIGUEATZONS FOR COMPLEXES NZTH A 
PARTIALLY PILLED M SHIU 


Complex 

Configuretlon 

Terms 

(l'2V) 

l»^2*2pSi* 

^8 


3«3d 

’d 


3P* 

^8, ®P, ^D 


3d* 

^8. ®P. S. ®P. 'g 


l»*2»*2p*3i*3p‘3d® 

'S. ®P. *D. ®P, 


3.* 3d^0 

's. ®P. S 


, 3.3p®3d^° 

'd. ®D 


! 3p^3d^° 

^8 

(I’aV)" 

la*2»*2p®3s3p 

'P. ^P 


3p3d 

*P. ®P. *D. ®D. 'P, ®P 


le*2i*2p‘3i*3p*3d® 

^P, ®P, *D, ®D, *F, ®P 


3. 3p»3d^® 

'P. ^ 







TABLE D-2 


wj In (1^3^). 

Is the 

term*dependent 

part of 

Wi - Wi + wj . 

(Wi)c f denotes 

is neglected. 

the value 

of obtained when configuration mixing 

Configuretlon 

Term 


(«i>c.f. 

V(«l>c.f. 

3s2 

's 

-.00357 

0 

-.00357 

3s3d 


-.00228 

-.00228 

- 



.01581 

.00684 

.00897 

3P^ 


-.00481 

-.00432 

-.00049 



-.00560 

.00432 

-.00992 


^8 

.01748 

.01727 

.00021 

3d2 


-.00564 

-.00564 

- 


'd 

.00295 

.00201 

.00094 



.00373 

.00322 

.00051 



.00616 

.00616 

- 



.02719 

.02382 

.00337 



TABLE D-3 


in (lV3^)° (See Table D-2) 


Conflg. 

Term 

«1 

(Wj) c.f. 

W^-(Wj)c.f. 

3r3p 


.01801 

.02116 

-.00315 



-.00753 

-.00705 

-.00048 

3p3d 


.02215 

.01900 

.00315 


3p° 

.00255 

.00207 

.00048 


'd° 

-.00869 

-.00869 

- 


V 

.00086 

.00086 

- 



.01934 

.01934 

- 



-.00859 

-.00859 

- 
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TABLE D-4 




In (1^3^^. 

(See Table D-2) 


Configuration Term 


<Wl’>cf 


3p^ 3dl0 

's 

.00262 

0. 

+.002620 

38 3p^ 3d® 


-.00228 

-.00228 

.... 



.01926 

.00684 

+.012424 

3*2 3p" 3dl0 


-.00406 

-.00432 

+.000258 


Id 

-.00763 

.00432 

-.011951 



.01625 

.01727 

-.001024 

38^ 3p6 3d8 


-.00564 

-.00564 

.... 



.00154 

.00201 

-.000472 



.00296 

.00322 

-.000258 


'g 

.00616 

.00616 

.... 


's 

.02226 

.02382 

-.001595 
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TABLE D-5 


In (1^2®3^*). (See Table D-2) 


Configuration Term 

“l* 



3s 3p5 3dl° 


.02281 

.02116 

.001649 



-.00678 

-.00705 

.000274 

38^ 3p^ 3d^ 

lp° 

.01735 

.01900 

-.001649 


3p° 

.00180 

.00207 

-.000274 



-.00869 

-.00869 

.... 



.00086 

.00086 

.... 



.01934 

.01934 

.... 



-.00859 

-.00859 

.... 


cf 



TABLE D-6 


THE UNITARY MATRICES j j^rF | tl>| 



3-2 


3p^ 

3d^ 

3.^ 

.99031 


.13789 

-.01670 

3p^ 

-.13490 


.98347 

.12079 

3d^ 

.03308 

is. 

-.11737 

.99254 


38 3d 


3p^ 

3d2 

3t3d 

.39843 


.43093 

-.01285 

3p^ 

-.43732 


.89702 

.06405 

3d^ 

.39642 

i. 

-.05193 

.99786 

raSEIH 


3p* 


3d2 

3p^ 


.99884 


.04807 

3d^ 


-.04807 


.99884 
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TABU D-7 


TW UNITMt MAniCI>||<^tr|c^||,<l’l*lV 



uyp 

3»3il 


3«3p 

SpM 

3«a» 

.Mia 

.1287S 


.99877 

.04956 

3»3d 

•.UB7S 

.MIM 

3pM 

-.04956 

.99877 


I 

I 


\ 


m 







TABLE D-8 


THE UHITARY MATRICES | |<trr|tf>j | 




3.*3pSd“ 

7?^ 

3p^3d^° 

.99491 

-.09904 

.01833 


.09991 

.99345 

.05552 

l.‘3,‘3d* 

-.01271 

.05707 

.99829 





6 9 

3*3p®3d* 

3.^3p^3al® 

3.Vjd* 

3»3|»‘3d® 

.831773 

.555116 

.000637 

3»*3p*3d^® 

-.554800 

.831337 

-.032791 

3.*3pV 

-.018733 

.026921 

.999462 



3.23p"3d'« 

3«^3p^3d* 

3»*3p^3d^® 

.999707 

-:024220 

3t^3p*3d^ 

.024220 

.999707 
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TABLE D-9 


THE UNITARY KATRICES | |<irr | tr>l | 


^ —. 



3.23p53d® 

3s3p^3d^° 

.99770 


-.06782 

3s^3p^3d^ 

_ 

.06762 


.99770 

3P° 


3s3p^3d^° 


38^3p^3d^ 

3s3p^3d^° 

.99960 


-.02818 

38^3p^3d^ 

.02818 


.99960 
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We call attention to the following features of the results. 


(i) In the complexes considered here, the effects of con¬ 
figuration mixing within a given complex are comparable in mag¬ 
nitude with the corresponding effects in complexes of the form 
(1^2^) studied previously. Thus, the mixing coefficients, which 
indicate the magnitude of departures from the central field 
approximation, are of the order of .1 and the changes in the 
intervals between adjacent terms are of the order of 10-30%. 

(ii) In a few instances the changes may be much greater 
than the values mentioned in the last paragraph. For example, 
the configurations 3s3d and 3p^ in the complex (1^2 3^) are 
very strongly mixed, the mixing coefficients being about .44. 

The same ia true of the corresponding pair of configurations 
in the complex (1^2^3 °), where the mixing coefficients are 
about .55. In these cases the resulting change in the value 
of W. is great enough to change the original order of the 
spectroscopic term (for large Z). Thus W^( D) actually lies 
below the value of W^(^P) in the configuration 3p^ of (1^2^3^) 
when allowance is made for configuration mixing within the 
complex. The change in W^ for this term thus exceeds 100% 

of the separation » Wj(^D) - Wj^(^P) as predicted by the 
Slater theory. 

(iii) It is necessary in general to allow for all effects 
of configuration mixing within a given complex, though in parti¬ 
cular cases soma ccMitributions may be negligible. In particular, 
.interactions connecting configurations that differ in the quantum 
numbers of a single electron cannot in general be neglected be- 
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aid* Intaractions that connect configurations that differ through 
the quantun nundDers of two electrons. The configurations 3s3d 
2 

and 3p discussed in the last paragraph Illustrate this point. 

(iv) The effects of configuration mixing are smaller in 
the complex (1^2^3^^) than in (1^2^3^). This beheirior is similar 
to that observed in complexes of the type (1^2^)t The effects 
of configuration mixing tend to diminish with increasing q. The 
reason is that the off-diagonal elements of electrostatic inter¬ 
action are essentially independent of the number of electrons in 
the system, whereas the differences between diagonal elements 
of the energy matrix increase with increasing H (total nusiber of 
electrons). Since the mixing coefficients are, in a first approxi¬ 
mation, equal to the ratios bet%ireen off-diagonal interactiMi 
elements and diagonal energy differences, they will tend to diminish 
with increasing N. 

P.l. The Smmi-Bapirical Method of Bacher and Goudsmit 

The method of Bacher and Goudsmit seeks to relate 

the term energies of an H-electron atom with the term energies 
of its ions, so that if the appropriate term energies of systMM 
(M,Z) with N < N are known (from either theory or experiment) 
one can calculate the term energies of the system (N,Z) to a 
certain approximation. 

The theoretical relations connecting the N-electron ener¬ 
gies with the 1-, 2-, ... , (N - 1)-electron energies are based 
on the ii^pendent-orbltal approximation. That is, each electron 
is assigned a definite pair of orbital quantum nusibers nl. We 
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now derive these relations« using a more explicit notation 
than that of Bacher and Goudsmlt. 

Let denote the nundser of Is orbitals, the number of 
2s orbitals, q^ the number of 2p orbitals, and so on up to q^, 
so that for an N~electron atom 

**1 **2 • • • + " • 

Let denote the set of quantum numbers that are required. In 
addition to the q^, to specify the state of an N-electron atom. 

Consider a given state An M-electron subsystem, 

consisting of the nucleus and M < H electrons, will not In 

1 M 1 

general be In a definite state (a ... a )r^. The a denote 
pairs of orbital quantum nusbers, not necessarily all distinct, 
selected from the set of N pairs that characterise the N- 

electron state under consideration. We shall use the abbrevl- 
1 M 

at Ion (a ... a ) « Thus we nay use the symbol a||r^ in 

place of (qj^ ... With every state we may associ¬ 
ate a definite probability * Here generally, we 

consider the set of all probabilities of the type > 

where A < B < N and the set Is contained In the set a^. 

These probabilities are not all Independent, since they evi¬ 
dently satisfy all relations of the form 

^p(r^irBmrglrc) - . (o-i) 

B 

(For sls^llclty we have omitted the labels a.) 

By means of what are now standard methods In the theory of 
complex spectra, one can evaluate any probability • 
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is the mean energy associated with the subconfiguration 
when the N-electron system is in state Note that the ener¬ 
gies do not coincide with the energies of the cor¬ 

responding M-electron ions of the N-electron system which would 
be denoted by , because the electrons in the sub¬ 

system move in the combined field of a nuclear charge Z and N-1 
orbital electrons, while the electrons in the ion move in the 
combined field of the same nuclear charge and M-1 orbital elec¬ 
trons. Thus the energy levels of the subsystem lie above the 
corresponding levels of the ion. However, in practical appli¬ 
cations of the method to be described, Bacher and Goudamit used 
empirical ionic energies as approximations to the corresponding 
subsystem energies. The latter cannot, of course, be obtained 
directly from experiment. We shall return to this question 
later. 


Following Bacher and Goudsmit, we write the energy of an 
M-electron subsystem as a weighted sum of contributions « 

K < M, where o„ stands for a set of K labels selected frcmi 
^ 1 M 

the set of M labels = « ... a . The energy contributions 

w are defined implicitly by the set of equations 
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(0-2) 



Th« Man energy contributions W are defined in analogy with 
the Man energies Wt 

•■k 

the set of all these equations, for all states otjil'u '*ith H < M, 
has a unique solution for the energy contributions w in terms 
of the energies Wt 


“m-1 “m-? 

- ... - (-I)" ^(1(0) 
a 


In this equation all energies pertain to subsystems of the 
N-electrMt system in the state proof of the fonmala 

is by MtheMtlcal induction. For M > 1, the formula reduces to 


W(a) « w(a) 


which coincides with equation (D-2) for M >■ 1. Suf^se now 
that the formula holds for M » 1, ... , n. We now show that 
it must then hold for M > n-i-1. Prom equation (D-4) we have. 


Vn 


(D-3) 


(D-4) 


(»-5) 
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X - I‘’K-/n-2lVn>*'<“n-2''n-2' 


Using the probablllty-ccsnpoBitlon rule (li-l)t, we car. put the 
right side of this relation into the form 




n+l-k^“n+l^n+l^ * 


The number of individual contributions to the kth summand is 
• This number is also given by 

k k-t-1 

j-O j-0 

Using the addition theorem for binomial coefficients, we obtain 


j-O 


Sc - 

which completes the proof. 


Bacher and Goudsmit assume that, for a given stats 
the energy contributions 


]... . 
(D-6) 

(D-7) 


(D-8) 


(D-9) 

( - 10 ) 
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diminish rapidly with k. Consider the expansion 


W(a^r^) - Y»<“) + Y * ••• * I '"‘“klVN’ * 

“ “2 “k 

* * “<Vh> ' 

which is a special case of equation (D-2). By omitting the 
last term on the right we obtain an approximate formula for 
the energy of the N-electron system in terms of quantities 
pertaining to k-electron subsystems with k ^ N. From equa¬ 
tion (D-4) it follows that the approximation * 0 is 

equivalent to the approximation 

"<Vh> ■ ■ 

”n-1 "n-2 “ 

Ulus the energy of the M-electron system is expressed as a 
linear cosdoination of energies of k-electron subsystems with 
k < N. 


More generally, we obtain the kth approximation in the 
sense of Bacher and Gioudsroit when we retain the first k-fl 
terns on the right side of equation ('0-11). I i particular, 
the first approximation (which is exact for 2-electron sys¬ 
tems) is defined by the equation 


^W(02|a^y - (M-2) 


(/O-ll) 


( 0 - 12 ) 


('0-13) 
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(k) 

With the help of equation (D-4) we can express W , for 
any k < Nj as a linear combination of energies pertaining 
to (k+l)-electron subsystems. 

In practical applications of the Bacher-Goudsmit approxi¬ 
mation method one uses experimental values of the ionic ener¬ 
gies place of the subsystem energies . 

In order for this to be a good approximation, the electrons 
that are present in the M-electron system but not in the k- 
electron subsystem must make only a minor contribution to the 
electrostatic field seen by the remaining electrons. Since 
the removal of "core electrons'* (electrons in closed shells) 
does substantially alter the effective electrostatic field of 
the atom, one should always regard the core electrons and the 
field they produce as fixed. That is, one should reckon all 
energies from the core energy as a fixed zero-point. The 
formulae already developed apply without change to these 
relative energies, since they do not depend on any assump¬ 
tions regarding the nuclear contribution to the electrostatic 
field except the assumption of radial symnetry. Since the 
core, by definition, is composed of coeqplete groups of elec¬ 
trons, its c<mtributi(m to the field will have spherical 
symmetry and may therefore be lusqped with the nuclear con¬ 
tribution. 

D-l:_The Bacher-Goudsmit Method and Configuration Mixing 

The independeht-orbital approximation is just as basic 
in the theory of Bacher and Ooudsmit as in the theory*of 
Slater. The basic equations of the theory, (D-2) and (D-4), 


213 


presuppose the possibility of assigning every subsystem of an 
M-electrcm system to a definite electronic configuration, and 
the coefficients in these equations are to be evaluated by 
means of a theory based on the same approximation. Neverthe¬ 
less, the Bacher-Goudsmit approximation does, in a sense, 
allow for certain effects of configuraticm mixing, because 
it makes use of empirical ionic energies, since these include 
contributions from configuraticm mixing, so roust the "predicted" 
energies. Now, the approximation of greatest practical and 
theoretical interest in the Bacher-Goudsmit scheme is the first 
one, which yields the energies defined by equaticm (1>-13). 

The question therefore arisest ' Which effects of ccmfiguraticm 
mixing are included in W^^^ when the two-electron energies 
^(g 2^2^°N^N^ taken from experiment? 

Bacher and Goudsmit discuss this question briefly and 
schematically, using ordinary perturbation theory. They con¬ 
clude (incorrectly, as we shall see) that the first approxi¬ 
mation includes all effects of configuraticm mixing with 
doubly excited configurati«ns (i.e., configurations differ¬ 
ing from the original one in two pairs of one-electron quan¬ 
tum numbers). They also conclude (correctly) that the first 
apfsroximation does not include all effects of configuration 
mixing with singly excited configurations. Before turning 
to a detailed consideration of these questions, however, let 
us examine the problem from another angle. 

As has previously been pointed out (Layser t)M 

Slater theory predicts certain qualitative regularities in 
the structure of term spectra, which do not accord with the 
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observed regularities. For example^ the Slater theory pre¬ 
dicts that all configurations consisting of n or 6-n equiv¬ 
alent p electrons, together with any number of complete 
groups, have the same structure apart from a scale factor 
Specifically, if n = 2, 3, or 4, the two independent inter¬ 
vals between the spectroscopic terms with L » 0, 1, 2 are 
predicted to satisfy the relation SD/DP = 3/2. The Z- 
dependent approximation scheme predicts different quali¬ 
tative and quantitative regularities. It predicts, for 
example, that the interval ratio SD/DP approaches a con¬ 
stant value with increasing Z in a given isoelectronic 
sequence, and that this constant (which can be calculated 
from the theory) depends both on the principal quantum 
numbers of the p electrons and on the complete groups pres¬ 
ent in the electronic configuration How does the Bacher- 
Goudsmit theory stand on this matter? 

2 2 2 

Consider a simple example, the configuration Is 2s 2p 
in the carbon sequence. We regard the K-electrons as com¬ 
posing a fixed core, and proceed to calculate the term 
energies relative to this core in the approximation (0-13). 

We obtain 

( 1 ) 3 2 2 -- 2 

W'^'(2s^2p SL) « W(2s ) + 4W{2s2p) + W(2p SL) - 4W(2s) - 4W(2p), (D-14 

where we have omitted unnecessary labels on the right side 
Mote that the mean two-electron energy W(2s2p) does not depend 
on the quantum numbers SL, because the two 2s electrons form 

a complete group. From this equation we see at once that the 
2 2 

configuration 2s 2p is predicted to have the same structure 
as the configuration 2p^. Actually, the interval ratio SD/DP 
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is substantially less than 3/2 in all systems belonging to 

2 2 

the beryllium sequence (i.e., in all configurations Is 2p ). 
The difference between the two limiting ratios is due to 
the mixing of configurations belonging to the same complex; 
both limiting ratios are correctly predicted by the Z-depen- 
dent theory. We conclude, therefore, that the Bacher*>6oudsmit 
method does not in general allow correctly for this type of 
configuration mixing — the most important type according to 
the Z-dependent theory — in the first approximation. 

Now, in the example just given the “perturbing" con¬ 
figurations that need to be considered if iche limiting inter¬ 
val ratios are to be correctly predicted are 2p (interacting 

2 2 2 2 
with 2s 2p ) and 2s (interacting with 2p ). These are 

“doubly excited" configurations in the sense explained above 

(though the configuration 2s^ should perhaps be described 

as doubly de-excited relative to 2p^). It follows that the 

first approximation in the Bacher-Goudsmit scheme does not, 

in fact, allow fully for the effects of mixing with doubly 

excited configurations. 

The argument given by Bacher and Goudsmit fails because 
it does not allow properly for the Pauli exclusion principle. 
Bacher and Goudsmit tacitly assume a strict correspondence 
between doubly excited N-electron configuraticms and doubly 
excited configurations in which the principal quantum num¬ 
bers of the excited electrons differ frcxs those of the orig¬ 
inal electrons, it does not, however, extend to doubly 
excited configurations with the same principal quantum num¬ 
bers as the original configuration. Thus, the de-excitation 
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of the 2p electrons in the configuration 2p^ to form the 
2 

configuration 2s is permitted by the exclusion principle; 

2 2 

but in the configuration 2s 2p , de-excitation of the 2p 

4 

electrons would yield the non-permitted configuration 2s . 

The predictions of the first-order Bacher-Goudsmit 
theory differ from those of the first-order Z-dependent 
theory in other important respects. For example, we have 
seen (see subsection 2) that the effects of mixing with 
singly excited configurations may be very large. We have 
also seen that in some instances the effects of configura¬ 
tion mixing are too large to be adequately described by 
second-order perturbation theory Finally, we have seen 
that the effects of configuration mixing are systematically 
greater in almost-eropty shells than in the corresponding 
almost-fullshells. All these features of the Z-dependent 
theory are at variance with the first-order Bacher-Goudsmit 
theory. 

We conclude that the Bacher-Goudsmit theory does not 
provide a satisfactory freunework for dealing with the most 
important effects of configuration mixing, even semi-em-. 
pirically. 

D .5. The Modified Slater Theory 

The theory of Bacher and Goudsroit is held by Racah 
to provide a theoretical justification for the 
L(L + 1) correction in the following sense. By means of 
equation (D-13) one can expreat the energies of an H-elec- 
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tron system in terms of the energies of two-electron sys¬ 
tems. Let us admit, for the sake of the argument, that the 
dominant effects of configuration mixing for the N-electron 
system are included in the energies calculated with equation 
(D-13) if one uses the true values of the energies ^ 
in evaluating the right side. Instead of using experimental 
values for the last-named energies, one could use an interpo¬ 
lation formula containing a number of adjustable parameters 
Suppose we are considering configurations of the type d” 

The configuration d^ has five terms, so that if one uses an 
interpolation formula with five adjustable parameters, the 
values of the parameters will be uniquely determined by 
experiment. Equation (D-13) then enables one to express 
the energy of any term belonging to the configuration d*' 
in terms of these five parameters* According to Slater's 
theory.^ the mean energy of all the terms in a given con¬ 
figuration labelled by given values of S and L is a linear 
combination of certain radial integrals, the Slater P and G 
integrals. For configurations d" there are three such 
integrals; P®, p^ and P^. However, even if one regards 
these as adjustable parameters, one cannot obtain a very 
good representation of the experimental data. The addition 
of two empirical correction terms, aL(L 1) 4- PQ, the secMid 


*In practice the parameters are not regarded as having 
fixed values, corresponding to the configuratim d^, but 
as functions of the effective central field. Onm can 
evaluate the parameters for n > 2 by equating the pre¬ 
dicted energies to their experimental counterparts and 
solving the resulting equations by the method of least- 
squares . 
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of which was suggested by Racah^^**^^ makes it possible to represent the 
2 

term energies in the configuration d exactly and to represent the term 
energies in other configurations that have been examined with fairly good 
accuracy. 

From the discussion in the preceding subsections, however, we know 
that equation (D-13) does not provide a satisfactory framework for dealing 
with the dominant effects of configuration mixing. It is true that the 
empirical description of spectra in the iron group based on the Slater- 
Racah interpolation formulae has been remarkably successful. However, 
guided by the considerations of the present paper, one can easily choose 
examples of configurations for which the method will give poor results. 

A simple example is the configuration ls^2s^2p^3p^. 
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APPENDIX E 
D. Layzer 

RELATIVISTIC EFFECTS IN HEAVY ATOMS 

E.1 Introduction 

In heavy, as in light, atoms the chief groupings of energy levels 
have their source In the nonrelatlvlstlc (electrostatic) Interactions 
between the electrons. These groupings are characterized by sets of 
prlnlcprl quantum numbers ( Adjacent groups are separated by 

Intervals of the order of a^mc^, where is the effective nuclear 
charge for the outermost electrons. 

The structure of a given group, however, has an entirely different 
character In heavy atoms .'n<I in light atoms. In the latter, each group 
(or complex) splits up Into spectroscopic terms oSL consisting of levels, 

dlstl igulshed by the quantum number J. The spacing between terms Is of 
2 2 

the order of Z^a me , while the spacing between levels in a given term 

4 4 2 2 3 

Is of the order of Z^a me • In light atoms, by definition, a Z^ « 1, 

so that the fine structure Is clearly differentiated from the coarse 

structure. The fine structure originates In relativistic Interactions 

between the orbital electrons and the nucleus, t* c course .strnv turt in 

nonrelatlvlstlc Interactions between the electrons. 

2 3 

In heavy atoms, on the other hand, the coupling parameter a Z^ Is 
not small. Consequently the relativistic Interactions cannot adequately 
be treated by first-order perturbation theory, as they can be In light 
atoms. The structure for the major groupings (complexes) Is correspondli^ly 
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complicated, in general, there Is no clear-cut distinction between 
coarse structure and fine structure; the two kinds of interactions con¬ 
tribute comparably to the spacing between levels in a given complex. 

3 2 

vnien Z^a » 1> the relativistic contributions become dominant; 
the complex breaks up into sub complexes , characterized by quantum 
numbers (n^j^). Within a given sub complex the spacing of levels is 
due primarily to the nonrelativistic (electrostatic) interaction. 

This state of affairs is characteristic of x-ray spectra. Here 
may be comparable to Z itself because the missing electrons belong to 
inner shells. 

The conventional method for describing the structure of energy 

2 3 

spectra in the intermediate-coupling regime (q 1) is set forth in 
(E-1) 

Theory of Atomic Spectra. It is based on the central-fieId ap¬ 

proximation and allows for the leading relativistic contribution to 
the energy - the spin-orbit interaction. This is evaluated for electrons 
moving, not in the Coulondi field due to the nuclear charge Z, but in an 
effective central field due to the nuclear charge plus the electrons. 

The conventional method has several drawbacks. 

It is better suited to a seml-eaq>irical description than a 
purely theoretical one. This is because the radial integrals ( () . 
that determine the swgnitude of the spin-orbit interaction cannot be 
accurately evaluated by purely theoretical means. They are best 
evaluated by treating them as adjustable parameters, to be assigned 
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valuaa that will make the theoretical formulae fit the data as well 
aa poaalhle. Thla la nearly alwaya the adopted In practice. 

The apln-orblt Interaction appeara aa a aeparate term In the 
Haailltonlan for a many-electron atom only In the Pauli (low Z) ap- 
proKlaMtloa. Whan Z la not aauill, It la by no meana clear that the 
apln-orblt Interaction adequately approxlmatea to the relatlvlatlc 
part of the Hamiltonian. 

An alternative deacrlptlon of relatlvlatlc effecta In the 

atructure and apectra of many-electron atoma la provided by recent work 
(E-2) 

of Layaer and Bahcall .■ Thla la an extanalon of the nonrelatlvlatlc 
Z-dependent theory^®"^^ The principal featurea of the Z-dependent 
relatlvlatlc theory are aa followa: 

(a) The Hamiltonian eaiployed la ahown to contain all known 

relatlvlatlc Interactlona up to (but not Including) thoae of order 
4 2 2 

a Z me . The largeat known Interactlona not Included In the Haodltonlan 
are thoae reaponalble for the lead> ahlft. 

(b) The apln-orblt Interaction la allowed for exactly, In 
the aeroth approximation - Juat aa In the Dirac theory of the hydrogen 
atom. 

(c) The theory enablea one to evaluate all contrlbutlona of 
order Z*^ coa^ared with the apln-orblt Interaction. Thla la not poa- 
alble with the conventional theory - and herein Ilea the crucial dif¬ 
ference between the old and the new methoda. It la the poaalblllty 

of preclaely evaluating the next-to-the-leadlng relatlvlatlc contrlbutlona 
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that renders Che new nechod quenCiCeClve where Che old one wee narely 
seal-quenCicaCIve. 

AlChough Che nexC-Co-Che leading relecivlaclc conCribuCions can be 

evelueCed In principle, Che ecCuel celculacions ere long end difficulc. 

2 

Ivan Cor such eoaperecively siaple elecCronic configurscions as Is 2s 
end ls^3s, considered by Leyser end Behcell^^ Che eveluacion of Chase 
Cema required a long celculecion. For heavy aCoaie Che corresponding 
ealeuleCiona are ec presenc prohibicively difficulc; IC is Cherefare of 
inceresc Co inquire whecher chare is en approxiaecion cnac will enable 
one Co esciaece Che conCribuCions in quesCion wich fair accuracy. Vm 
deal wich Chia quaacion in Saccion B.2. In SecCion E.2, we consider a 
relacad problaa (as ic Curns ouc Co be)-Che ralaciviscic generalisacion 
of Che nonralaciviscic screening paraaecera defined by Layser^^*^^'. 




2 4 2 

The coafficienc of a Z ac in cha eKprasaion for Che energy of an 

acoaic scace, according co Che Cheory of Layser and Bahcall^''^''^^ 

2 3 

is a funccion of cha coupling paraascer X ■ a Z . In cha oxcreaM rela- 
Civlscic case, chia coefficienc, which we denoce by a^, is given by a 
direcc generalisaclon of Soaaaarfeld'a fonaila: 



In Che opposite llaic (X 0) e^ does £££ i^xoach the value that would 
be predicted on the basis of the central-field epproxiaationi to Ataln 
the correct Halting valua^ one aust allow for the effects of configuration 
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mixing. (This can be done most simply by using, in place of the central* 
field wave function \|r, a finite linear combination of wave functions 
]^a^\|r^, the constants a^ being given by the screening theory; see 
Appendix C. 

In the general case a^(X) can be calculated, without excessive 
labor, to any desired accuracy by methods described in the paper by 
Layser and Bahcall cited earlier. We may, therefore, consider it to 
be known. 

4 3 2 

It is the coefficient a^(X) of the term a Z me that presents the 
real problem, and which we would like to approximate. 

Now in the standard theory of one>electron spectra (including x*ray 
spectra), the relativistic contribution to the energy is usually 
written in the form 

a^(Z - 0)^ a^mc^ + 0(aV) (E-2) 

where a^ is given by Sommerfeld 's formula, 

and a is a suitably chosen constant. Evidently 



Thus a knowledge of o is equivalent to a knowledge of a^. 

Now soaw, though not all, of the contributions to a^ can be inter* 
preted physically in terms of screening. For example, tlM spin-other- 
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orbit interactions, which contribute to a^, tend to reduce the magnitude 
of the (true) spin-orbit interaction in much the same way as the inter- 
electronic electrostatic Interactions reduce the magnitude of the electro¬ 
static inter .ctions between the nucleus and the electrons. On the other 
hand, such two-electron interactions as the spin-spin interaction have 
no one-electron counterpart. The interpretation of a as a screening 
parameter is therefore an oversimplification. Nevertheless, it may 
furnish a useful means of estimating the value of a^ through Equation (r>3). 

The simplest assumption one can make about a is to equate it to S, 
the screening parameter that figures in the screening theory. A com¬ 
parison of emprlcal values of o for hydrogenlc spectra with theoretical 
values of S shows that this a a fair approximation, though not as good 
as one might wish. 

Fortunately, the fractional errors due to inaccuracies in the 
estimated value of o tend to become progressively smaller with increasing 
Z because the quartic term becomes progressively larger in comparison 
with the cubic term. Thus, the screening approximation to a^ should 
be considerably better for Uranium, say than for Lithium. 

A semi-empirical study to improve the estimate of a might proceed 
along the following lines: 

The two-electron screening parameter S(u/3) [representing the 
screening of an a electron by a ^ electron] Is given In a first ap¬ 
proximation by the formula^^'^^ 
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s(a/p) 


I j Pp (rS dr^ 


(E-4) 


This can be written as 


S(a/^) -j wCr^) P ^r^ < r^ dr^ (E-5) 

i.e. as a certain weighted average of the probability P that 

a 3 electron will be found closer to the nucleus than r , the distance 

a 

of an a electron. The weight w Is given by 

This weighting factor is appropriate for electrostatic screening, but 
not necessarily for other kinds of screening. A more general law than 
(E-6) is: 

"k ■ 4 <’a> 'a'" 

It would be interesting to test whether, by an appropriate choice of k, 
one could arrive at a formula for o in better agreement with experiment. 
[The connection between the total screening constants and the two- 
electron screening parameters is discussed in the paper by Layzer^^ 
cited earlier. ] In the meantime, the old definition (S ■ o) can be 
used with a reasonable expectation of fair accuracy. 
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E.3 Relativistic Screening Parameters 


The screening parameters that we have dealt with so far are defined 
by a nonrelatlvlstlc theory. The question naturally arises as to what 
constitutes the proper relativistic generalization of the screening 
parameters. 

To begin with, the relativistic screening parameters must depend 
on the coupling parameter X as well as on Z. They must, in fact, have 
the form 

S(X.Z) - S^(X) + S_j(X)z‘^ + 0(Z‘^). (E-8) 

In the limit X -* 0, S(X,Z) -♦ the nonrelavistic screening parameter S(0,Z). 
The coefficients S^(X) are all bounded from above and below. In the 
relativistic limit X -• oo, 

S(X,Z) - S^(oo) (E-9) 

since the remaining terms on the right side of (a) vanish in this 
limit. 

It is possible to develop a detailed theory of the X-dependent 
screening parameters along the lines of the variational nonrelatlvlstlc 
theory. The main result is that the contributions to the screening 
paranwters from electrons in closed shells are independent of X« 

Now the X-dependence of the screening parameters could be significant 
only in heavy atoms. Since most of the electrons in a heavy atom are 
norswlly found in closed shells, the amplitude of variation of the 
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screening parameters with X must be small, so that the nonrelativlstlc 
approximation to the screening parameters should be adequate for most 
practical purposes. 
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APPENDIX F 

D. Layzer and A. Naqvl 

STATISTICAL DESCRIPTION OF COMPLEX SPECTRA 

F. 1. Mean Energies of Configurations and Complexes 

In any partial description of complex energy spectra the 
mean energies of configurations and complexes play a central 
role. Even if no other information about the energy spectrum 
is available one can still obtain from the mean energies, with 
the help of Saha's equation, useful estimates for the concen¬ 
trations of the various atomic and ionic species in a mixture 
of gases. Again, in constructing radial wavefunctions for 
use in the evaluation of absolute transition .probabilities — 
either by the screening-approximation or by a more refined 
method — one may not wish to use a different set of wave- 
functions for each distinct level. In a first approximation 
one could then use a single set of wavefunctions correspond¬ 
ing to the mean energy. 

If we neglect entirely the effects of configuration mix¬ 
ing, we can write the mean interaction energy of any configu¬ 
ration as a sum of two-electron contributions (nl,n'l'). If 
the atom has N electrons the mean interaction energy consists 

of (j) » N(N - l)/2 such contributions. For example, the 

^ 2 2 5 

mean energy of the configuration Is 2s 2p is 

(Isis) -t- (2s2s) 4 10(2p2p) + 4(ls2s) + 10(ls2p) -i- 10(2s2p) . (F-1) 

Here N ■ 9 and the number of two-electron interactions is 
9 X 8/2 » 36. 

The two-electron interactions can be expressed in terms of 
the Slater F and G Integrals. The appropriate formulae are 
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given by Layzer and Slater but are In^llclt 

in, earlier work by Condon and Shortley and Shortley. 


To calculate the total mean energy of a configuration, 
including the one-electron part, we nay conveniently eBg>loy 
the Z-dependent screening theory. The mean energy is then 
given by the formula 

“ - - I • 

Here is the number of electrons of the type nl in the 
configuration under consideration and 

2 , » Z - s , . 
nl nl 

The screening parameters are defined in terms of the 
t%#o-electron interaction energies: 

■nl ■ - *nl,nT>*<"*l"'‘'> ' 

n*l' 


n^d(nl,n’I' )/Sz^^ (nl i<n'l') 

s(nl|n'l') • 

in^d{nl,nl)/^z^^ (nl -n'l') 

Ihe two-electron screening paraateters depend on Z« approach¬ 
ing CMistancy as z •» «. 


(F-2) 


(F-3) 


(F-4) 


(F-5) 
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The preceding discussion can easily be extended to com¬ 
plexes. The advantage of doing this is that the resulting 
formula for the mean energy of a complex will have the cor¬ 
rect linear dependence on Z. On the other hand, if one is 
dealing with a complex that contains a large number of con¬ 
figurations, the extra work required to allow properly for 
all the interactions may offset the resultant gain in accuracy. 

F. 2. The Coefficients for Complex Configurations 

Instead of working with the energies themselves, one may 
wish to consider the statistical properties of the energy co¬ 
efficients W., defined as the coefficient of Z in an expansion 
of the energy in powers of Z In this subsection we discuss 
briefly some relevant properties of the W^, before returning, 
in the following subsection, to the main stream of our discus¬ 
sion. The considerations of the subsequent subsections will 
apply equally to the energies themselves and to the coeffici¬ 
ents W^. 

The coefficient is composed of two parts. The two- 
electron part can be calculated from the two-electron 
interactions and is the same for all levels belonging to a 
given configuration. The remainder of W^, which we denote 
by W^, results partly frtxn LS-dependent interactions within 
the configuration and partly from interactions connecting 

terms belonging to different configurations in the same com- 

(F~l) 

plex. It has been shown (Layzer ) that in atoms with a 
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partly-filled L shell these two contributions to are often 
comparable in magnitude Some recent, as yet unpublished, cal¬ 
culations show that the same is true of atoms with partly-filled 
M shells Hence the effects of configuration mixing within a 
complex cannot in general be neglected beside the LS-dependent 
contributions from interactions within a given configuration. 

As a preliminary to calculations of thai allow properly 
for the effects of configuration mixing, and as a temporary sub¬ 
stitute for such calculations, we may conveniently evaluate the 
contributions to from interactions within the configuration. 
This is by no means a trivial matter for complex configurations, 
since such configurations may contain several distinct spectro¬ 
scopic terras characterized by given values of S and L. In the 
representations normally used the energy operator is diagonal 
with respect to S and L but not with respect to the quantum 
numbers that serve to distinguish spectroscopic terms with the 
same S and L Hence the energy matrix, instead of being dia¬ 
gonal. will in general consist of muTti-diiu'jnsional blocks. 

This IS the form in which energy calculations have usually been 
left in the past. In order to complete the diagonalization 
one needs to know the ratios of the .various Slater integrals 
that figure in the energy matrix elements In the standard 
theory of complex spectra these ratios are regarded either as 
parameters to be determined semi-empirically. or as theoretical 
quantities whose evaluation requires the solution of a set of 
Hartree-Fock equations On the other hand in the Z-dependent 
scheme one evaluates the Slater integrals using straight hydro- 
genic wavefunctions the contributions to fzam interactions 
within a given configuration have definite numerical values 
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which can be calculated once and for all, as can the elements 
of the transformation matrices that diagonalize the blocks re¬ 
ferring to spectroscopic terms of given S and L. The trans¬ 
formation to a representation in which has no nondiagonal 
elements connecting terms in the same configuration is the 
first step toward the ultimate objective of a representation 
in which the energy matrix has no nondiagonal elements connec¬ 
ting terms belonging to the same complex. 

P. 3. The Variance of the Energy 

The variance of the energy (in a given configuration) is 
defined by 

var W = <(W - W)^> g 0 ^ . (F-6) 

Moszkowski has written down formulae that give for 

an arbitrarily complex configuration in terms of the o^s for 
two-electron configurations. The latter can easily be evalu¬ 
ated by direct calculation. Proofs of Moszkowski's formulae 
have not as yet appeared; but for the purposes of this report 
we shall assume that the formulae are correct — or at least 
that 0 ^ can be found for any configuration of interest. 

The parameter measures, roughly, the width of a con¬ 
figuration. Knowing a^, one can estimate (from Boltzmann's 
formula) the variation in the population of levels over a 
given ccmfiguration. More generally, one can estimate the 
error that results, in any particular calculation, from 
assigning the same energy to all levels in a given configu¬ 
ration. From one can also derive a rough estimate for 
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the average spacing between levels in a configuration. Thus 
if one assumes that, say, 50% of the levels lie in the range 
W ± 0 ^, one can reasonably expect the resulting estimate of 
the level spacing to be accurate to within a factor of 2. 

F. 4. The Distribution Function of the Levels in a Configuration 

Knowing o we know the width of the probability distri¬ 
bution, but not its shape. The question now arises: What 
is the best assumption as to the shape of the distribution 
that is consistent with our limited knowledge of the energy 
spectrum? 

In addition to W and o we know from direct calculations 
w 

that the distribution of energy levels in a configuration is 
strongly skewed. This fact has suggested (unpublished Convalr 
report, 1961) that the distribution can be represented by an 
exponential function of the energy. Once the exponential form 
has been chosen, the two parameters on which such a distribu¬ 
tion depends are uniquely determined by the values of W and o^. 

In fact, the distribution is given by 

N exp [-(w - W)/o^l (W > W - o^) 

P(W) •= (P-7) 

0 (W < W - o^) 

Thus the energy spectrum is supposed to have a sharp cut-off 
0 ^ units below the mean energy W. The probability distribution 
achieves its maximum value at the cut-off and then falls off 
exponentially. 
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Before this law, or any such law, can be compared with experiment 
we need to have some reasonable and consistent way of passing from the 
theoretical or experimental energy spectrum, which is discrete, to a con¬ 
tinuous probability distribution. One such way is the following. We 
assign to each level the appropriate statistical weight (2S + 1). Then 
we divide the energy range into a convenient number of equal intervals, 

bounded on the left by the lowest energy in the configuration, and 
determine the total statistical weight of all levels falling in each 
energy interval. We then plot a histogram of the running means of the 
statistical weights against the energy. This corresponds to smearing 
out the discrete energy levels, which occur in bundles, into a somewhat 
continuous distribution of energy. As long as the intervals used to 

smooth the data are reasonable small compared with the total energy 
range, the exact manner of smoothing is of no importance. One should 
bear in mind that, as discussed in Subsection F.2, the distribution of 
energy levels as calculated with neglect of configuration mixing is at 
best semi-quantitative, so that much of the detail present in an un¬ 
smoothed plot of the distribution of energy levels is spurious anyway. 

Figures F-1 thru F-6 illustrate the smoothed distributions of 

energy levels in configurations consisting of equivalent d and f elec- 

(F-8) 

trons. Experimental data from the compilation by Miss Moore , are 
taken for the d-electron configurations and theoretical energies are 

taken for the f-electron configuration from the paper by Elliot, Judd 
(F-9) 

and Runclman. If in the experimental data all theoretical 

terms are not Identified, the missing terms are given in each 
figure. These plots exemplify the skewed character of energy dis¬ 
tribution, but Che skewness appears to be much less marked 
than that associated with the exponential 
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Figure F-1. Smoothed Distribution of Experimental Energy 
Levels in V II (Configuration 3d^) 
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Figure F-2. Smoothed Distribution of Experimental Energy 
Levels in Cr II (Configuration 3d^) 
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Figure F-3. &KX>thed Dletrlbucion of Experimental Energy 
Levels in Nb II (Configuration 4d^). 
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Level in the Configuration 



distribution. In every case the main peak of the distribution 
lies significantly to the right of the lower cut-off. 

Another defect of the exponential distribution consists in 
its not having the lower cut-off in the right place. The lower 
cut-off has in fact been arbitrarily placed units below the 
mean energy W. This assignment is usually in error by a quan¬ 
tity not much smaller than o itself. Thus the main features 
* w 

of the distribution are not as well represented by the expon¬ 
ential law as one would like. 

We have sought to improve the statistical description in 
two ways. Firsts we have added a third parameter, which can 
in every case be calculated exactly with comparatively little 
labor. Secondly, we have chosen the distribution function in 
accordance with an objective principle — the principle of 
maximum entropy inference — which ensures that the probability 
distribution incorporates the maximum possible uncertainty 
consistent with having the correct mean, variance, and lower 
cut-off. These two improvements are discussed in detail in the 
following subsections. 

F. 5. The Lower Cut-off 

The lower cut-off is simply the energy of the ground 
state of the configuration. Now, according to a well-known 
rule of Hund, the lowest spectroscopic term in a given con¬ 
figuration is the one with the highest value of S; or. If two 
or more terms have the same value of S, the one of these with 
the highest value of L. One of the levels belonging to this 
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term, namely, the level with » S and » l, is charac¬ 
terized not only by definite values of M_, S, M., and the one- 

i ^ S 

electron azimuthal quantum numbers I , but also by definite 
values of the one-electron magnetic quantum numbers m^. Hence 
the corresponding wavefunction is a Slater determinant, and 
the corresponding energy can be obtained by a straightforward 
calculation. As an example we calculate the energy of the 
ground state of the configuration f^ . (See Section F-9) 


F. 6. Maximum Entropy Inference 

Consider the problem of determining the probability dis¬ 
tribution p(x) associated with a given random variable X. 

Suppose that the possible values of X are known to lie in the 
range (a,b) and that the first k moments of X are also knownt 

b 

4^ - xS(x)dx (i - 1, ... ,k) . (P-8) 

a 

p(x) must satisfy the further conditions 
b 

P(x) >0, f p(x)dx - 1 . (F-9) 

a 

According to the principle of maximum entropy inference (Jaynes 
the best choice for the function p(x) is the one that 
maximises the entropy 

b 

H - - J p(x) in p(x)dx , (P-10) 
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subject to the k + 1 auxiliary conditions expressed by the 
above integral equations. Using the method of Lagrange multi¬ 
pliers^ one easily finds that p(x) must have the form 


2 k 

p(x) = N exp [a^x + a^x + ... + a^x | . 

The constants a, ... a, , N, can be found from the k + 1 inte- 
1 k 

gral conditions, 

The entropy of a probability distribution is a measure of 

(F-7) 

the uncertainty inherent in the distribution. Shannon 
showed that no other function can consistently bear this inter¬ 
pretation. Thus the principle of maximum entropy inference 
may fairly be described as the correct method for choosing a 

probability distribution when partial information about the 

(F-6) 

distribution exists. As Jaynes has pointed out, it is 

a quantitative embodiment of Laplace's principle of insuffi¬ 
cient reason. 

F. 7. Application to the Di st ribution of Energy Levels in a Con ¬ 
figuration 

We now have three parameters at our disposal — the mean 
energy W, the standard deviation of the energy o^, and the 
energy of the ground state, which we denote by W^. Let 

X = (W - Wq)/€ , 


(F-11) 


(F-12) 


I 


where a is an as yet unspecified unit of energy. The variable 
X is a dimensionless measure of the energy relative to the 
ground state. 


Applying the principle of maximum entropy inference, we 
find that the probability distribution has ihc form 

p = N exp - b)“i) (x > 0) . 

The energy unit e has been chosen so as to make the coefficient 
2 

of X in the exponential -1. By integrating the last equation 
over the range of x we obtain a formula for the normalization 
constant N; 

N = 27r"^[l + erf(b) ]"^ 
where erf(x) is the error function 

, ^ 2 

erf(x) = 2 Tr “2 e"^ dy . 

0 

Thus p(x) depends on the single parameter b. 

Let 


. - (w - Wq)A , a « aye ; 

u and a are the mean energy and the standard deviation in 
dimensionless units. 

The integral conditions which ensure that the probability 
distribution has the required mean and variance yield two inde¬ 
pendent relations fr<»n which b and e can be evaluated in terms 
of u and 0 . The following form of these rSlatioi)s is conveni¬ 
ent for calculationt 

U « b + ip(0) , 

(0A) • » 11 - np(0) ]/2u^ . 


(F-13) 


(F-14) 

(F-15) 


(F-l6ab) 


(F-17) 

(P-18) 
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By substituting the first equation into the second, we obtain a trans¬ 
cendental equation for b in terms of the ratio (a/(i) . The numarical 
solution of this equation for a range of values of a/|i can be obtained 
easily. Table F-1 gives a/ii for various values of b, and tlM sasM is 
plotted in Figure F-7. Having obtained the value of b from this 
table, one can evaluate the energy unit e from equations (F-17) and 
(F-16a). 

The agreement between the predicted and observed distributions 
seems to be satisfactory. If future studies should disclose a 
systematic discrepancy between the predicted and observed distribu¬ 
tion, then the very existence of such discrepancies would represent 
fresh information about the probability distribution and could be 
Incorporated into the conditions of the problem. Moreover, it is 
likely that any such systematic discrepancy would have an assignable 
physical cause. Thus the principle of maximum entropy Inference 
appears to be a powerful method for Isolating those properties of com¬ 
plex spectra that are susceptible to interpretation in reasonably 
siflq>le terms. 

f iL. tgmiptiqn gt. tht 99iblLMU9ft Swttrf 

Consider a transition array, defined as the aggregate of lines 
resulting from all possible transitions connecting two 
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b 

S. 

-- 

Jti- 

00 

.75122 

C.l 

.70515 

0.25 

.63495 

0.5 

.59937 

0.75 

.51272 

1.0 

.53512 

1.25 

.48832 

1.5 

.43851 

1.75 

.39107 

2.0 

.34742 

2.5 

.28235 

3.0 

.23567 

3.5 

.20203 

4.0 

.17677 
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subject to the k + 1 auxiliary conditions expressed by the 
above integral equations. Using the method of Lagrange multi¬ 
pliers, one easily finds that p(x) must have the form 

2 )c 

p(x) = N exp [a^^x + a^x + ... + ' 1 • 

The constants a, ... a, , N, can be found from the k -f 1 inte- 
1 k 

gral conditions. 

The entropy of a probability distribution is a measure of 

(F-7 ) 

the uncertainty in^ierent in the distribution. Shannon 
showed that no other function can consistently bear this inter¬ 
pretation. Thus the principle of maximum entropy inference 
may fairly be described as the correct method for choosing a 
probability distribution when partial information about the 
distribution exists. As Jaynes ^ has pointed out, it is 

a quantitative embodiment of Laplace's principle of insuffi-- 
dent reason. 

F. 7. Application to the nistribution of Energy Levels in a Con ¬ 
figuration 

We now have three parameters at our disposal — the mean 
energy W, the standard deviation of the energy o^, and the 
energy of the ground state, which we denote by W^. Let 

X » (W - W^)/e , 

where e is an as yet unspecified unit of energy. The variable 
X is a dimensionless measure of the energy relative to the 
ground state. 


m 


(F-ll) 


(F-12) 
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Figure F-7. Solution to Equations F-17 and F-18. 



given configurations. What can we say about the probability 
distribution of line frequencies in such an array, given the 
probability distributions for the energy levels of the two 
configurations? 


Let denote the energy of the ith level in the configu¬ 
ration a. The frequency corresponding to a transition ai bj 
is given by Bohr's frequency condition, 

hv(ai,bj) . 

We assume for the sake of definiteness that configuration a 
lies above configuration b, so that the frequencies defined 
by this equation are always positive. 

If departures from LS coupling are not too great, the 
combination spectrum [the spectrum of frequencies v(al,bj)] 
will have a clumpy structure. Transitions between levels 
belonging to the same pair of spectroscopic terms give rise 
to multiplets . and in LS coupling "the separation bet%ifeen 
multiplets is large compared with the separation between 
lines in a multiplet. 

The distribution of lines in a multiplet and the dis-. 
tribution of multiplets in a transition array present distint 
problems. In a first approximation we may assume that the 
two distributions are stochastically independent — that is, 

2Sl 


(F-19) 



that the form of the probability distribution for lines in a 
multiplet does not depend on the position of the multiplet in 
its transition array. Similar remarks apply to the case of 
pure jj coupling. 

On the other handj if the coupling scheme is not close to 
either LS or jj coupling, neither the individual energy spectra 
nor the combination spectrum will exhibit any marked fine struc¬ 
ture. In this respect the case of intermediate.coupling is 
easier to deal with than either extreme LS coupling or extreme 
jj coupling. 

There is a second important respect in which the stochastic 
description of the combination spectrum is simpler for inter¬ 
mediate coupling than for LS coupling. In LS coupling, transi¬ 
tions connecting levels belonging to terms of different multi¬ 
plicity (different S) are strictly forbidden. This well-known 
selection' rule substantially reduces the number of lines in a 
transition array that occur with nonvanishing intensity. More¬ 
over, it forces one to deal separately with the various classes 
of transitions corresponding to the various possible values of S. 
In jj coupling transitions connecting terms differing by more 
than one quantum number j are forbidden. In intermediate coup¬ 
ling, on the other hand, intersystem transitions are not only 
allowed, but often have high probability, since intermediate 
coupling prevails in heavy atoms, we may reasonably assume that 
the choice of an initial level and the choice of a final level 
are ato<diastically independent events. The remainder of this 
discussicHi proceeds from this assumption. 
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If we regard the quantltiea v, as random variables 

[whose possible values are the numbers v(al,bj)re¬ 
spectively then the aaaunqptlon just stated means that the 
random variables and are stochastically Independent^ so 
that the random variable v Is the difference between two stoch¬ 
astically Independent random variables. Olils fact enables us 
to write down the probability distribution of v In terms of the 
probability distributions of and According to a well- 

known theorem In probability theory, the first distribution Is 
just the convolution of the probability distributions of and 
W^. Let Pj^Cx) denote the probability distribution of and p^ 
denote the probability distribution of Then the probability 

distribution p^^ of 


(F -20) 


Is given by 


Pab^y^ “ - x)dx . (F-21) 

From this equation, or, more directly, from the definition of 
stochastic Independence, one can show that the expectation and 
the variance of are given by 


(f-22) 


Var W . « Var W, + Var W. . 


(F-23) 


As an Illustration of these rules we consider the example 
of two normal distributions. By an appropriate choice of sero- 
points we can always ensure that « 0, so that we 
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may put 


p (x) » (27rcr^)"^ exp [-x^Sa^] 


p^(x) = (27ra^)"2 

Since the normal probability distribution ia syimnetric^ the 
probability distribution Pj^ of coincides with p^ in the 
present example; but this is not true in general. A short 
calculation, in which one performs the integration by ccxn- 
pleting the square in the exponent, gives 

p ^(x) = (2-rro\ exp [-xV2o^. ] , 


4 




the probability distribution of the difference is also 
normal and its variance is the sum of the variances of and 
W^, This is the well-known "reproductive property" of the 
normal distribution. 

We assume that the probability distributions p^, have 
the form derived earlier in the present report. It is conveni¬ 
ent to introduce the function 

1 (X > 0) 

€(x) - 

0 (X < 0) 

With the help of this functi<Mi we can write the probability 
distribution for W. in the form 


P_(x) « <P(x;a .p )€(x) , 


(F-24) 

(F-25) 


(F-26) 

(F-27) 


(P-28) 

(F-29) 
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where 


<p(x;o^P^) = N^exp [-(x - ^^)^/2a^] . ( 

The definition of p^ is exactly analogous. Note that 

5^(x) - p^(-x) ( 

We therefore obtain for the probability distribution of W^^^ 


**ab^y^ " f " y;Oj^>Pjj) e(x)€(x - y)dx 


From the definition of the function €(x) it follows at once, that 


€(X)€(X - y) 


e(x - y) (y > 0) 


€(x) (y < 0) 


**ab^^^ “ J (y < 0) ; 

0 

00 

Pab<y) = J' ‘P(xfO^,p^)<p(x - y;Ojj,Pjj)dx 

y 

00 

J <P(y + d* (y > 0) . 


The integrations are straightforward (ms the Matheaaticsl Mote at 
the end of this (^psndlx). The final result is: 


F-30) 

F-31) 

F-32) 

F-33) 

F-34a) 

F-34b) 
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Pab<y> “ 


( y < 0) 


/P J-, ^ + y <r 

1 - erf ( :-sa,'„i* r 


(F-35a) 


- P.,(y) 


where P^^(y) . ^ 


b^a ^ y ^a \ 

(STzT^ / 

< 

Va®^b "^2 + <rj^2) ^ (y > 0) 

(P-35b) 


and erf(x) has been defined in equation (F*15) 


F.9. Ground State Bnarav for a Configuration of Kauivalent Electrons 

He consider a configuration of equivalent electrons Hm energy 
of the coaq>letely filled shells is consaon to all the tersv of the con¬ 
figuration, It can be calculated, if desired, but is of no conse¬ 

quence for our purpose. FOrtheraore, the average energy of such a 

(F-2) 

configuration can be calculated in a straightforward manner (Slater ). 

W,v (i^) (P-36) 

2i 

and (i.i) P‘‘ (ni, nX) (F-37) 

k 

where Che coefficiente f for conflgurecioiis of p, d» f end g electrone 

V 

are given in Table J-1 and F are the Slater P integrals defined in 
Equations (B-26) and<B-27). 
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Us shall describe below a method of calculating the energy of the 
ground term. 

It Is well known that a subshell i can contain at most 2 ( 2 £ + 1) electrons 
and Is therefore exactly half filled when q ■ 2i + 1. It Is also well 
known that the term energies (with respect to the average energy) of a 
configuration with n electrons are the same as those of a configuration 
with [2(2i + 1) " q) electrons. It Is, therefore, necessary only to 
consider configurations which are half or less than half filled, l.e. 

< 2f + 1. 

The ground term of any configuration has the largest value of S 
for that configuration. If several terms have the same largest value of 
S, say It Is the one with the largest value of L, say 

asKmg them. Each term Is highly degenerate, but one of the SL 
states belonging to the term (S^^, Is such that 

\ • ‘•«x. 

The sero*order state, which belongs uniquely to the above set of values 
of the quantum nusdiers can be found easily. Mg > 8 ^^ arises idien all 
spins are parallel and In this case of all n electrons Is +1 (.r -a) 
Therefore, 

Mg » - i q (f-38) 

and the multiplicity of the ground term Is (2 + 1) or (q 1). 
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Tlie largest value of consistent with the above requirement on the 
and Is obtained by a slnq>le application of the exclusion principle. 
The values of the q equivalent electrons must be 

= /, i - I, .. (i - q f 1) (F-39) 

For this case 

\ ‘ f «' • -t * 1) 

The corresponding zero order state Is designated In the usual way by 
(f)' (i-l)"^^. (/ - q + 1)^ 


where the quantity In parenthesis Is the value of of Individual 
electrons and the sign Indicates that m^ -■ 1-1/2. Tlie ground terms of 
a few configuration of type i" are given In Table P-2. For a shell wlilch 
Is exactly half filled (q - IL * 1), It Is readily seen from Equation 

(F-5) that s 0 and the term Is ^ 

TABLE F-2 

Ground Terms for Configurations of Equivalent Electrons 
Configuration Multiplicity L Ground Term 

1= q + 1] I- J (2i-q+l)J 


P 

d 

f 

g 


2 

3 

7 

5 


3 

4 
8 
6 


3 

0 

10 


P 

^’f 

«S 

6k, 
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The energy of the ground term in the absence of external fields is 
the same as that of any one of the zero order states of this term, say 

of the zero order state mentioned above. This energy can be found by 
(P-3) 

standard methods (Condon and Shortley ), and Is given by, 

W -V y I a** (i m^, i iS^) - b^ (i m“, / mj)! p'^ (ni, ni) 

(P-41) 

where the first auiaaatlon la over all two electrons combinations and 
the second susmatlon Is over k which takes on all even values from zero 
to 2t. The coefficients a^ nnd b*^ are defined by Equations 8^ (14) and 

8^ (6) of Condon and Shortleyand are tabulated by them In Tables 
6 6 

1 and 2 for s, p, d and f electrons. 

^ As an example we calculate below the energy of the ground term 

I of configuration f^. The zero order state of lowest energy is 

3"^ 2 * 1^ O'^ -1^ -2^ -s’*" 

and uniquely belongs to the term. Table P-3 gives the appropriate 
k k 

e and b for f electrons. From the sums of these quentities, given in 
the last row of this table, it readily follows that the energy of the 
ground term, ^8, of the configuration f^ is 

1 
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TABLE F-3 


®i» '* £-electrons 

Mote: e^ « 1 for ell esses 


b® » 0 If ^ ai^ 


Electrons 

a p 

m 

X 1089 

m 

: b2 

1 

X 225 

b^ 

X 1089 

b® 

X 7361.64 

3"^ 2'" 

0 

-21 

-6 

25 

30 

7 

1^ 

-15 

3 

15 

10 

54 

28 

0+ 

-20 

18 

-20 

0 

63 

84 

-1^ 

-15 

3 

15 

0 

42 

210 

-2^ 

0 

-21 * 

-6 

0 

0 

462 

-3^ 

25 

9 

1 

0 

0 

924 

2"- 1+ 

■I 


-90 

15 

32 

105 

0^ 

mm 


120 

20 

3 

224 

-1^ 

n 

H 

-90 

0 

14 

378 

-2* 


:^9 

36 

0 

70 

504 

-3^ 

Bl 


-6 

0 

0 

462 

0* 

12 

6 

-300 

2 

15 

350 

-1* 

9 

1 

225 

24 

40 

420 

-2* 

0 

-7 

-90 

0 

14 

378 

-3* 

-15 

3 

15 

0 

42 

210 

0"^ -l^ 

12 

6 

-300 

2 

15 

350 

-2^ 

0 

-42 

120 

20 

3 

224 

- 3 * 

-20 

18 

•20 

0 

63 

84 

- l ^ - 2 * 

0 

-7 

-90 

15 

32 

105 

- 3 ^ 

-15 

3 

15 

10 

54 

28 


mm 

-21 

-6 

25 

30 

7 

n 

m 

-77 

-462 

168 

• 

616 

5544 


2b& 

























^•10. Mathaaatlcal Itof 


From Bquationa (F-30) and F-34a) wa obtain, for y < 0, 


00 

r 

(X - p^)^ (X - y - p^)^ ■ 

’\hj “P 

2 a ^ 2 (j. ^ 

O 

La b J 

00 

r ^ 

■“b y •"p 

-(bk^ + 2qx + r) dx 

«• J 

o 



dx (P-43a) 


(P-43b) 


<T 2 + a ^ 

, -A_b— 


(P-44a) 


. p, * (y » 

' -a^V 

* (y * 


(P-44b) 


(P-44c) 


1h« InCagral on tha right hand alda of Equation (P*43b) can bn avaluatad. 
by maana of formula (bb), TabU 314, paga 65 of tha conpilation of 
daflnlta Intograla by GrSbnar and Hofraitar^^'^^^ Applying this formula 
wa obtain Equation (P-35a). 


Por y > 0, «M combino Equations (P*30) and (P-34b) and gat 



f r <2 - (-y ♦ @.) 

P.t(y) -W J «p--I'*- 

o L a 


(2 - 


(P-45) 


which Is the same as Equation (P-42a) If we Interchange the subscripts 
a and b and replace y by (>y). Performing the integration as before 
we obtain Equation (P-3Sb). 
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APPENDIX G 

H. E. Stubbs and G. A. Victor 

THE EFFECT OF LINE WIDTH, STRENGTH, AND SHAPE 
ON RADIATIVE TRANSFER IN AN EXPANDING GAS 

G.l Introduction 

The details of line width, shape, and strength have a very important 
role in the calculation of radiative transfer phenomena in non-expanding 
gas clouds. This dependence on details of line shape introduces very 
serious complications in most astrophysical applications^^ since 

reliable information on line shapes is available only for simple gases 
as hydrogen and heliumThe object of this appendix is to investi¬ 
gate a simple auxiel of a rapidly expanding gas which results in expressions 
for line absorption which are largely independent of detail line shape. 

The origin of the simolifieation in this model is found in the fact that 
the Doppler shift due to the velocity of expansion becomes the dominant 
factor, and that the details of line shape due to collision and thermal 
Doppler broadening become unimportant. 

TItis appendix will be divided into three sections. Part G.2 will 
be devoted to some simple estimations of line widths due to the thermal 
Doppler broadening, and to compare it with the Doppler shift ^e to the 
expansion velocity. Part G.3 will show that resonance absorption is 
independent of line shape. Part G.4 will present a simple calculation 
for an ei^nding hydr^an gas for the transmission of the 82, 260 em~^ 
line. 
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G.2 Doppler Shift Du« to V<loclty of Expanalon 


I 


We shell construct the following ideellsed model for our expending 
cloud: 


Am nusitber deneity of hydrogen etoms is given hy: 

p (r.t) - 1.5 X 10^/t^ perticles/cm^ r < R^(t) 

p (r.t) - 0 r > Il^(t) 

R^(t) - V^t 

- 1 X 10*^ c (6-1) 

t is the time in seconds 
c is the speed of light > 3 x lO^^cm/sec. 


the reletive velocity between two perticles is elong their line of 
centers end is given by: 

(C-2) 

Hhere r^2 perticles seperetion: 

The reletivistic Doppler formule is: 


X • X^p cos Q (6-3) 

Since Che velocity is elong the line of centers we cen set cos 6 ■ 1. 
Since « C we heve the reletion > .(t2. therefore: 

Sua ■ *0*^0 

The fmnmtle for the thenMl Doppler width is 
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aXj - X 


(G-5) 



for Hydrogen we cen construct the following table: 


3/2 kT 


V W 

100 ev 

3.2 X 10*\ 

1.6 X 10** 

10 ev 


5.0 X lO"* 

1 ev 

3.2 X 10’\ 

o 

1.6 X lO'* 

.1 ev 

lo’h 

5 X 10*^ 


o 



Froa Table G*>1 we see that the Doppler width due to therael aotion is 
always aaich ssutller than the sHUciwai Doppler shift due to the velocity 
of expansion. The results of the rest of this appendix depend on this 
fact. 


G.3 Independence of Resonance Absorption on Line Shape 

Let a photon of frequency w be enitted at r ■ 0. Let the 
probability of emissimi of a photon of energy between hu end If (w -f du) 
per unit tiaa per particle be giyen by P(wX The probability that a given 
aton will absorb a photon of energy Vu is given by, 

p(«) - k P(w) (G-6) 

from the principle of detail balance. 
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The probability that a photon of frequency (t> will tranaverae a dig' 


tance r before being absorbed is given by: 

J A(|,w)d| 

I (u,r) ■ e ° 


(G-7) 


with 


A(6,u) - p(|) p (u*[r,u ]) 


(G-8) 


where p(t) is the density of particles which can abaorb this resonance 
radiation. I(Utr) can therefore be written as: 


f 


P(t) P (4i‘ 


(G-9) 


l(u,r) ■ e 
^ere we have used the Doppler foraula: 

a (G*>10) 

and evaluated this at a particular tisM t^. 

We intergrate this over all frequencies, weighted with P(u) to got the 
probability that any photon will go a distance r. 


T(r) 


•f 


+00 

du P(u) e 




idiich will be re-written as 


I<r) 


p® .k r dio 

ml dw P(w) e 


<l) ^1) 


(G-11) 


( 0 - 12 ) 
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If p(|) It a constant out to a distance r, which is consistent with our 
cloud nodal, we have: 


T(r) 


■ / du P(w) ^ 

JmCD 


u*Au(r) 

Hi)di 


(G-13) 


where: K* * kp Aa)(r) • 


Finally 


«(!- 


+ f K‘P(u')du* 

P(w)du e 

•U) 


Using the result of part G.2 we can calculate this for r ■ oo 
*00 

-+00 

T(a»)-/ P(«)dwe*' 


•£ 


J^’ju »(«')<•“' 


r“ 

if we let T) ■/ P(a)')du' 
then dr) ■ F(u)du 


(C-14) 


(6-15) 


T(oo) 




(6-16) 


fron Equation (6-16) we see that for an infinite cloud 
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resonance absorption is Independent of line shape. For e finite cloud, 
it is only the photons thetsre eaitted near the edge that will show resonance 
absorption that is dependent on the form of P(u). The ratio of the thick¬ 
ness of this edge to the cloud radius is essentially Ve there¬ 

fore eoncludle that resonance absorption is controlled priasurily by line 
strength, for this simple model of an expanding cloud. The object of 
this appendix is to investigate a siaqple model of a rapidly expanding gas 
^ich results in expressions for line absorption which are largely 
independent of detail line shape. 

Arguments similar to thoae .given above for resonance absorption can 
be given for line absorption due to two close, but non-overlapping lines. 

It is easy to see that an absorption line centered at can only absorb 
radiation from lines centered at such that: 

^ -““o 

It is again obvious that this absorption is independent of line 
shape except for lines where the equality holds in Equation (G-17). Aie 
effect of a finite cloud on the i^sorption of non-overlapping line radio- . 
tion does little but introduce some simple geometrical factors. 

G.4 Iftatmlision of. Jyd,roien 82,260 cm"^ ,fcing 

^ line of interest is the 82,260 cm*^ line which is nsitted by a 
transition from the first excited state to the ground state. We will 
assume that the gas has a temperature history as shoun in Figure 6-1. 

Xn addition, we assume that the gas is at thermal equilibrium at all times. 
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TEMPERATURE (tv) 


010 * 40-1910 



Figure G-1. Teaperaturc During Expansion. 
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At equilibrium, the population of etatee is given by the Saha equation 


■ number density of protmia 

> number density of electrons 

■ -onisation potential ( r degree) 

Q ■ principle quantum number 

a nuid)er density of atoms in state n 

'he number density of particles which can absorb the 82,260 cm'^ line is 
given by 

P(t,T) - (t)f(t,T) (C-19) 

NjCT) 

Istimations of f(t,T) from Equation(G.xg)lead to the results shown in the 
following table. 


Time 

fefltt^lys] 

R [cm] 

*t 


r 

00 

io“^ 

1.5 X 10*^ 

2 X 10* 

7 X 10*f 

5.7 X 10^ 

2.5 X 10'* 

10'* 

1.5 X lo'* 

2 X 10* 

7 X 10*f 

1.08X 10'^ 

1.3 X lO'* 

10'^ 

1.5 X 10^* 

2 X 10^ 

7 X 10*f 

.26 

1.4 X 10'^ 

10-’ 

I.JKIO" 

2 X 10* 

7 X lO^f 

1 

1.4 X 10'* 

10'* 

1.5 X 10® 

2 X 10^ 

1 7 X lO'^f 

1 





From Equation (G>16) vk hava: 


T^(co).i (l-.-'O 

t 

trtiere the subscript indicates time dependence, 
ing ei9ression: 


(G-20) 

is given by the follow 


. .0 . P (tjT) A t (G. 

^ 8jt V 8« 

o 

A is the Sinstein coefficient S x 10^ sec*^. The results for 
at various times is shown in Table G*2. FiguresG>2 end G-3 show the 
variation of T(ao) with time for this simple siodel. 


272 



FRACTION TRANSmtTTCD TO OOHT) 



Figure G-2. TrensaiMion in Expending Mndiua. 


I 
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